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Abstract. We rigorously justify in 3D the main asymptotic mod- 
els used in coastal oceanography, including: shallow-water equations, 
Boussinesq systems, Kadomtsev-Petviashvili (KP) approximation, 
Green-Naghdi equations, Serre approximation and full-dispersion 
model. We first introduce a "variable" nondimensionalized version of 
the water-waves equations which vary from shallow to deep water, 
and which involves four dimensionless parameters. Using a nonlocal 
energy adapted to the equations, we can prove a well-posedness theo- 
rem, uniformly with respect to all the parameters. Its validity ranges 
therefore from shallow to deep-water, from small to large surface and 
bottom variations, and from fully to weakly transverse waves. 
The physical regimes corresponding to the aforementioned models 
can therefore be studied as particular cases; it turns out that the 
existence time and the energy bounds given by the theorem are al- 
ways those needed to justify the asymptotic models. We can therefore 
derive and justify them in a systematic way. 



1. Introduction 

1.1. General setting 



The motion of a perfect, incompressible and irrotational fluid un- 
der the influence of gravity is described by the free surface Euler (or 
water-waves) equations. Their complexity led physicists and mathe- 
maticians to derive simpler sets of equations likely to describe the dy- 
namics of the water-waves equations in some specific physical regimes. 
In fact, many of the most famous equations of mathematical physics 
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were historically obtained as formal asymptotic limits of the water- 
waves equations: the shallow-water equations, the Korteweg-de Vries 
(KdV) and Kadomtsev-Petviashvili (KP) equations, the Boussinesq 
systems, etc. Each of these asymptotic limits corresponds to a very 
specific physical regime whose range of validity is determined in terms 
of the characteristics of the flow (amplitude, wavelength, anisotropy, 
bottom topography, depth, ...). 

The derivation of these models goes back to the XlXth century, but 
the rigorous analysis of their relevance as approximate models for the 
water-waves equations only began three decades ago with the works 
of Ovsjannikov [30llll], Craig [l3], and Kano and Nishida [271[28l|26] 
who first addressed the problem of justifying the formal asymptotics. 
For all the different asymptotic models, the problem can be formu- 
lated as follows: 1) do the water-waves equations have a solution on 
the time scale relevant for the asymptotic model and 2) does this 
model furnish a good approximation of the solution? Answering the 
first question requires a large-time existence theorem for the water- 
waves equations, while the second one requires a rigorous derivation 
of the asymptotic models and a precise control of the approximation 
error. 

Following the pioneer works for one-dimensional surfaces {IDH) 
of Ovsjannikov [H] and Nalimov [38] (see also Yosihara [50ll51j ). 
Craig [T3], and Kano and Nishida [27] provided the first justifica- 
tion of the KdV and IDH Boussinesq and shallow water approxi- 
mations. However, the comprehension of the well-posedness theory 
for the water-waves equations hindered the perspective of justifying 
the other asymptotic regimes until the breakthroughs of S. Wu ( [48j 
and [49j respectively for the IDH and 2DH case, in infinite depth, 
and without restrictive assumptions). Since then, the literature on 
free surface Euler equations has been very active: the case of finite 
depth was proved in [29j, and in the related case of the study of the 
free surface of a liquid in vacuum with zero gravity, Lindblad |34ll35j . 
More recently Coutand and Shkoller ^12j and Shatah and Zeng |45j 
managed to remove the irrotationality condition and/or took into ac- 
count surface tension effects (see also [3] for IDH water-waves with 
surface tension). 

In order to review the existing results of the rigorous justification 
of asymptotic models for water-waves, it is suitable to classify the 
different physical regimes using two dimensionless numbers: the am- 
plitude parameter e and the shallowness parameter ^ (defined below 
in 

— Shallow- water, large amplitude (/x <C 1, e ~ 1). Formally, this 
regime leads at first order to the well-known "shallow-water equa- 
tions" (or Saint- Venant) and at second order to the so-called 
"Green-Naghdi" model, often used in coastal oceanography be- 
cause it takes into account the dispersive effects neglected by 
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the shallow- water equations. The first rigorous justification of the 
shallow- water model goes back to Ovsjannikov |40ll41j and Kano 
and Nishida [27j who proved the convergence of the solutions of 
the shallow-water equations to solutions of the water-waves equa- 
tions as /i — > in IDH, and under some restrictive assumptions 
(small and analytic data). More recently, Y. A. Li [33j removed 
these assumptions and rigorously justified the shallow- water and 
Green-Naghdi equations, in IDH for flat bottoms. Finally, the 
first and so far only rigorous work on a 2DH asymptotic model is 
due to a very recent work by T. Iguchi [23] in which he justified 
the 2DH shallow-water equations, also allowing non-flat bottoms, 
but under a restrictive zero mass assumption on the velocity. 

— Shallow water, medium amplitude (/i ^ 1, e ~ ^'T'-)- This regime 
leads to the so-called Serre equations, which are quite similar to 
the aforementioned Green-Naghdi equations and are also often 
used in coastal oceanography. To our knowledge, no rigorous result 
exists on that model. 

— Shallow water, small amplitude (^u ^ 1, e ~ //). This regime (also 
called long-waves regime) leads to many mathematically interest- 
ing models due to the balance of nonlinear and dispersive effects: 

— Boussinesq systems: since the first derivation by Boussinesq, 
many formally equivalent models (also named after Boussi- 
nesq) have been derived. W. Craig [T3] and Kano and Nishida 
[28] were the first to give a full justification of these models, 
in IDH (and for flat bottoms and small data). Note, however, 
that the convergence result given in [28j is given on a time scale 
too short to capture the nonlinear and dispersive effects specific 
to the Boussinesq systems; in [13], the correct large time exis- 
tence (and convergence) results for the water-waves equations 
are given. The proof, of such a large time well-posedness result 
for the water-waves equations, is the most delicate point in the 
justification process. Furthermore, it is the last step needed 
to fully justify the Boussinesq systems in 2DH, owing to [5] 
(flat bottoms) and [9] (general bottom topography), where the 
convergence property is proved assuming that the large-time 
well-posedness theorem holds. 

— Uncoupled models: at first order, the Boussinesq systems re- 
duce to a simple wave equation and, in IDH, the motion of 
the free surface can be described as the sum of two uncoupled 
counter-propagating waves, slightly modulated by a Korteweg- 
de Vries (KdV) equation. In 2DH and for weakly transverse 
waves, a similar phenomenon occurs, but with the Kadomtsev- 
Petviashvili (KP) equation replacing the KdV equation. Many 
papers addressed the problem of validating the KdV model 
(e.g. [13.28.43.l5ll47t[23]) and its justification is now complete. 
For the KP model, a first attempt was done in [26j, under re- 
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strictive assumption (small and analytic data), but as in [27j . 
the time scale considered is unfortunately too small for the rel- 
evant dynamics. A series of works then proved the KP limit 
for simplified systems and toy models [20ll3ll32], while a differ- 
ent approach was used in [32] where the KP limit is proved 
for the full water-waves equations, assuming a large-time well- 
posedness theorem and a specific control of the solutions. 
— Deep-water, small steepness > 1, s^/JJ■ <C 1). This regime leads 
to the full-dispersion (or Matsuno) equations; to our knowledge, 
no rigorous result exists on this point. 

Instead of developing an existence/convergence theory for each 
physical scaling, we hereby propose a global method which allows 
one to justify all the asymptotics mentioned above at once. In order 
to do that, we nondimensionalize the water-waves equations, and keep 
track of the five physical quantities which characterize the dynamics 
of the water-waves: amplitude, depth, wavelength in the longitudinal 
direction, wavelength in the transverse direction and amplitude of 
the bottom variations. 

Our main theorem gives an estimate of the existence time of the so- 
lution of the water-waves equations which is uniform with respect 
to all these parameters. In order to prove this theorem, we intro- 
duce an energy which involves the aforementioned parameters and 
use it to construct our solution by an iterative scheme. Moreover, 
this energy provides some bounds on the solutions which appear to 
be exactly those needed in the justification of the asymptotics regimes 
mentioned above. 



1.2. Presentation of the results 

Parameterizing the free surface hy z = C(t, X) (with X = (x, y) G M?) 
and the bottom hy z = —d + h{X) (with c? > constant), one can 
use the incompressibility and irrotationality conditions to write the 
water-waves equations under Bernouilli's formulation, in terms of a 
velocity potential (j) (i.e., the velocity field is given by v = Vx,20): 



where V = {dx,dy)'^ and dn4> is the outward normal derivative at the 
boundary of the fiuid domain. 

The qualitative study of the water-waves equations is made easier 
by the introduction of dimensionless variables and unknowns. This 
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requires the introduction of various orders of magnitude linked to the 
physical regime under consideration. More precisely, let us introduce 
the following quantities: a is the order of amplitude of the waves; A 
is the wave-length of the waves in the x direction; A/7 is the wave- 
length of the waves in the y direction; B is the order of amplitude of 
the variations of the bottom topography. 
We also introduce the following dimensionless parameters 

-^=e, - = /?; (1.2) 

the parameter e is often called nonlinearity parameter, while /i is the 
shallowness parameter. In total generality, one has 

(e, fi, 7, (3) G (0, 1] X (0, 00) X (0, 1] X [0, 1] (1.3) 

(the conditions e £ (0, 1] and (3 G [0, 1] mean that the the surface and 
bottom variations are at most of the order of depth — /3 = corre- 
sponding to flat bottoms — and the condition 7 € (0, 1] says that the 
X axis is chosen to be the longitudinal direction for weakly transverse 
waves). 

Zakharov [52] remarked that the system (jl.ip could be written in 
Hamiltonian form in terms of the free surface elevation and of the 
trace of the velocity potential at the surface -0 = 4^1^-^ ^^id Craig, 
Sulem and Sulem [T7] and Craig, Schanz and Sulem [16j used the fact 
that (jl.ip could be reduced to a system of two evolution equations on 
and ip; this formulation has commonly been used since then. The 
dimensionless form of this formulation involves the parameters intro- 
duced in ()1.2p . the transversity 7, and a parameter v = {1 + y/Ji)~^ 
whose presence is due to the fact that the nondimensionalization is 
not the same in deep and shallow water. It is derived in Appendix [XI 

9.^ + C + ^|V-^l 2(l + eV|V^CP) = 

(1.4) 

where V''' = {dx,"fdy)'^ and Q f^^-yle^, f3b] is the Dirichlet-Neumann 
operator defined by Qf,^y[eC, I3b]^p = {I + e'^\VC\^)'^/^dn'P\^^^^, with <P 
solving 

+ fid^<P + 7V9y<^ = 0, -1 + Pb<z< eC, , s 

In Section [21 we give some preliminary results which will be used 
throughout the paper: a few technical results (such as commutator 
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estimates) are given in §2.11 and elliptic boundary value problems di- 
rectly linked to ()1.5I) are studied in §2.2[ 

Section [3] is devoted to the study of various aspects of the Dirichlet- 
Neumann operator Q^^y[eC, (3b]. It is well-known that the Dirichlet- 
Neumann operator is a pseudo-differential operator of order one; in 
particular, it acts continuously on Sobolev spaces and its operator 
norm, commutators with derivatives, etc., have been extensively stud- 
ied. The task here is more delicate because of the presence of four 
parameters (e, //,7,/3) in the operator G^^^[e(, (3b]. Indeed, some of 
the classical estimates on the Dirichlet-Neumann are not uniform 
with respect to the parameters and must be modified. But the main 
difficulty is due to the fact that the energy introduced in this paper 
is not of Sobolev type; namely, it is given by 



(1.6) 

where := (—9^ — 7^5^)^/^. For high frequencies, this energy is 

equivalent to the x i?*"'"^/^-norm specific to the non-strictly hy- 
perbolic nature of the water-waves equations (see [T4J for a detailed 
comment on this point), but the equivalence is not uniform with re- 
spect to the parameters, and the H'' x estimates of [48 1 149 1 
29t[3] for instance, are useless for our purposes here. We thus have 
to work with estimates in | • Ij^^-type norms and the classical re- 
sults on Sobolev estimates of pseudo differential operators cannot be 
used. Consequently, we must rely on the structural properties of the 
water-waves equations much more heavily than in the previous works 
quoted above. 

Fundamental properties of the DN operators are given in §3.H while 
commutator estimates and further properties are investigated in §3.21 
and §3.31 We then give asymptotic expansions of Q ^^j[e( , (3b]i(! in 
terms of the parameters in ^3.41 

Using the results of the previous sections, we study the Cauchy prob- 
lem associated to the linearization of (jl.4p in Section HI the main 
energy estimate is given in Proposition 14.11 

The full nonlinear equations are addressed in Section [5] and our main 
result is stated in Theorem 15.11 it gives a "large-time" (of order 
0{£/iy)) existence result for the water-waves equations (jl.4p and a 
bound on its energy (defined in (|1.6p ). The most important point 
is that this result is uniform with respect to all the parameters 
(e, /i, 7, (3) satisfying (|1.3p and such that the steepness £y/Jl and the 
ratio (3/e remain bounded. The theorem also requires a classical Tay- 
lor sign condition on the initial data; we give in Proposition [5T] very 
simple sufficient conditions (involving in particular the "anisotropic 
Hessian" of the bottom parameterization 6), which imply that the 
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Taylor sign condition is satisfied. 

Both Theorem 1 5 . 1 1 and Proposition 15.11 can be used for all the physical 
regimes given in the previous section, and the solution they provide 
exists over a time scale relevant with respect to the dynamics of the 
asymptotic models. We can therefore study the asymptotic limits, 
which is done in Section [6l It is convenient to use the classification 
introduced previously to present our results (we also refer to [31\ for 
an overview of the methods developed here): 

- Shallow- water, large amplitude (/i <C 1, e ~ 1). We justify in ^6.1.11 
the shallow-water equations without the restrictive assumptions of 
|24j and previous works. For the Green-Naghdi model, we extend 
in §6.1.2l the result of [33j to non-flat bottoms, and to two dimen- 
sional surfaces. 

- Shallow water, medium amplitude (/u <C 1, e ~ yT^)- We rigorously 
justify the Serre approximation over the relevant 0{\/ yfp) time 
scale in §6.1.21 

- Shallow water, small amplitude (/x <C 1, e ~ /u). 

— Boussinesq systems: In §6.21 we fully justify all the Boussinesq 
systems in the open case of two-dimensional surfaces (flat or 
non-flat bottoms). 

— Uncoupled models: We complete the full justification of the KP 
approximation in §6.31 

- Deep-water, small steepness (/i > 1, e^/Ji <C 1). We show in ^6.4.11 
that the solutions of the full-dispersion model converge to exact 
solutions of the water-waves equations as the steepness goes to 
zero and give accurate error estimates. 

We also give in ^6.4.21 an estimate on the precision of a model used 
for the numerical computation of the water-waves equations (see 
|15j for instance). 

1.3. Notations 

- We use the generic notation C(Ai, A2, . . . ) to denote a nondecreasing 
function of the parameters Ai, A2, . . • • 

- The notation a < b means that a < Cb, for some nonnegative 
constant C whose exact expression is of no importance (in particular, 
it is independent of the small parameters involved). 

- For all tempered distribution u £ we denote by u its Fourier 
transform. 

- Fourier multipliers: For all rapidly decaying u G ©(M'^) and all 
/ G C(]R^) with tempered growth, f{D) is the distribution defined 
by 

(this definition can be extended to wider spaces of functions). 

- We write (0 = (1 + and A = {D). 
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- For all 1 < p < oo, | • \p denotes the classical norm of L''(]R^) while 
II • lip stands for the canonical norm of L''(5), with S = M? x (—1,0). 

- For all s G M, H^{R^) is the classical Sobolev space defined as 

H'{R'^) = {ue e'im^), \u\h, := \A'u\2 < oo}. 

- For all s € M, II • ll^oo^i^s denotes the canonical norm of 
L°°([-l,0];i?^(M2)). 

- If B is a Banach space, then | • |bt stands for the canonical norm 

of L°°([0,T];S). 

- For all 7 > 0, we write V' = {dx,jdy)^, so that V' coincides with 
the usual gradient when 7=1. We also use the Fourier multiplier 
\D'^\ defined as 

as well as the anisotropic divergence operator 

div^, = (V^)^. 

- We denote by ^/<,7 (or simply !p when no confusion is possible) the 
Fourier multiplier of order 1/2 

- We write X = (x,y) and Vx,z = {dx,dy,dz)^; we also write 

- We use the condensed notation 

As = Bs + (C,),>, (1.9) 

to say that Ag = Bg if s < s and Ag = Bg + Cg if s > s. 

- By convention, we take 



l[pk = l and J2pk = 0. (1.10) 

k=l k=l 

- When the notation dnU^g^^ is used for boundary conditions of an 
elliptic equation of the form Vx,z ■ x,z'u = h \n some open set 
Q, it stands for the outward conormal derivative associated to this 
operator, namely, 

dnu\^n = ■ PV X,zU\^^, (1.11) 
n standing for the outward unit normal vector to dQ. 
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2. Preliminary results 

2.1. Commutator estimates and anisotropic Poisson regularization 

We recall first the tame product and Moser estimates in Sobolev 
spaces: if io > 1 and s > 0, then \/f e H^n H^o(j^2^yg g H'iR"^), 

Ifdln^ < \ f\H*o\g\H' + {\f\H'\9\H*o)s>to (2.1) 
and, for ah F E C"^(M";M"^) such that F(0) = 0, 

yu£H'{R^)'', F{u)£H'{R^r and \F{u)\hs < C{\u\oo)\u\hs ■ 

(2.2) 

In the next proposition, we give tame commutator estimates. 

Proposition 2.1 (Ths. 3 and 6 of [30j). Let to > 1 and -to < ^ < 

to+1. Then, for all s > 0, / G H^^+'^nH'+''{R'^) and u G H'+^''^\R'^), 

where we used the notation (|1.9p . 

One can deduce from the above proposition some commutator esti- 
mates useful in the present study. 

Corollary 2.1. Let to > 1, > and 7 G (0, 1]. Then: 

i. For all v G p H'+\R'^f and u G H'{R'^), one has 



[yl'',div^(v-)]n|2 < |v|^to+2|u|//s + (^\u\fjto+i\v\f^s+i] 



s>to+l' 



ii. For all < r < to + 1, / G L°°((-l, 0); n i7*o+^(IR2)) and 
u£ L^{{-l,0);H'+^-\R^)), 

Proof. For the first point, just remark that 

[yi",div^(v)]u = [yl^div^(v)]'u+ [yl^v] -v^u, 

and use Proposition 12.11 to obtain the result (recall that 7 < 1). 
For the second point of the corollary, remark that for all z G [—1,0], 



\A[A'J]uiz)L < \fiz)\Hto+i\uiz)\Hs + {\fiz)\Hs+i\u{z)\Hto] 



1>to 



as a consequence of Proposition 12.11 (with r = 1). The corollary then 
follows easily. □ 

Let us end this section with a result on anisotropic Poisson reg- 
ularization (when 7 = /i = 1, the result below is just the standard 
gain of half a derivative of the Poisson regularization). 
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Proposition 2.2. Let 76 {0,1], fi > and x be a smooth, compactly 
supported function and u G ©'(M^). Define also := x{\'l^z\D''\)u. 
For all s G M, if u G H''-'^/^{R'^), one has A^'u^ G L'^{S) and 

I 1 I „ t„ I 1 I 

Moreover, for all s G M, if u G H''+'^/^(M?), one has A''V''''^u^ G 
and 

^"^^1 ^.1 <||V1^V-V||2<C'J 



In the above estimates, ci, C2, c'^ and C2 are nonnegative constants 
which depend only on 

Proof. Write classically (with IC'I = VCi + 7^^!)' 

\\x{Vl^z\D^\)u\\l,= f f {0''x{VJiz\C\?\m\'dzdC 
7r2 J -I 



where F denotes a primitive of x^- The first estimate of the proposi- 
tion then follows from the observation that 

2 1 ^ m - 2 1 

where the constants depend only on x- 

For the second estimate of the proposition, remark that 

llV'^^'Tu^l^^ ^\\x{VJiz\D^\)mu\\^ + ^\\x'{z^\D^^^^^ 

and use the first part of the proposition. □ 



2.2. Elliptic estimates on a strip 

We recall that the velocity potential # within the fluid domain solves 
the boundary value elliptic problem 

r52# + /x52^ + 72^92^ = 0, -l + pb<z<eC, 

with (e,^,7,/3) G (0,1] x (0, 00) x (0,1] x [0,1]. 

Denote by 5 the flat strip 5 = M x (—1,0), and assume that the 

following assumption is satisfied: 

There exists ^0 > such that 1 + eC - pb > ho- (2.4) 
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Under this assumption, one can define a diffeomorphism S mapping 
S onto the fluid domain f2: 

S ^ Q 
^ ■ (X, z) ^ S{X, z) := {X, z + a{X, z)) , 

with 

a{X,z) = -(3zb{X)+e{z + l)CiX). (2.5) 



Remark 2.1. The mapping a used in (j2.5p to define the diffeomor- 
phism S is the most simple one can think of. If one wanted to have 
optimal estimates with respect to the fluid or bottom parameteriza- 
tion (but unfortunately not uniform with respect to the parameters), 
one should use instead regularizing diffeormorphisms as in Prop. 2.13 

of m. 



From Proposition 2.7 of [29], we know that the BVP (j2.3p is equiv- 
alent to the BVP (recall that we use the convention (II. lip for normal 
derivatives), 

/ Vx,^ • P[a]Vx,zcl> = 0, in 5, . . 

with (j) = <P o S and with the (2 -|- 1) x (2 -|- 1) matrix P[cr] given by 

Remark also that it follows from the expression of P[(t] that 

Vx,. • p[a]Vx,z = v^'T • (1 + gH)v^'^ 

where 

/ dzcr -^dxO- 

(2.7) 

Below we provide two important properties satisfied by Q[cr]. 

Proposition 2.3. Let to > 1, s > 0, and C,^ G n if''+^(M2) be 

such that (|2.4p is satisfied. Assume also that a is as defined in 
Then: 
i. One has 
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and, when a is also time dependent, 

1 

'ho 



ii. For all j > 1 and h G H*''+'^nH''+'^{M.'^y , and denoting by [a] h 
the j-th derivative of C ^ QW] i^i the direction h, one has 

||Q(^-)[ct] -hll^^^, < (£)^C(-3-,e^,||V'^'V||^^^.o) 
i j 

X ( I^Ih^+iH I^'Ih*o+i + ((1 + II V^'^alli^^HOn l^^l^^'o+i>,>,J • 
fc=l i^fc fc=l 

iii. The matrix 1 + Q[cr] is coercive in the sense that 

\0\^ <k[a]{l + Q[a])0-O, 

with 

1 / n2 

k[(j] := k^^^f[a] = 1 + ll^^o-Hoo + J^y^ + VmII^^^IIooJ • 

Proof. The first two points follow directly from the tame product and 
Moser's estimate ()2.ip and ()2.2p . and the explicit expression of Q[(7]. 
It is not difficult to see that {1 + Q[a])0 ■ = j:^\B0\^, where 

n+dzcr -y/Jld^a \ 
B=[ 1 + a^a -7^/if^5j,a . 
V 1 / 

The matrix B is invertible and its inverse is given by 

-, / 1 y/Jldxa \ 

1 + VoO 1 + 9^(7/ 

Remark now that owing to (j2.4p . the mapping o", as given by (j2.5p . 
satisfies (1 + dzcy)^"^ < that 



Vi + a^cj|B V3^K3 < Vi + ll^^f^lloo + ^(i + 



^0-||oo 



Since ^|]r3,_^r3 > I©!, the third claim of the proposition fol- 

lows. □ 
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Since the Dirichlet condition in (12. 6p can be 'lifted' in order to 
take homogeneous Dirichlet boundary condition, we are led to study 
the following class of elliptic BVPs: 

VM,7.(i + Q[^])VM.7n = V^'^-g, in 5, 
-"1,^0 = 0' 9nn|^_, = -e^-g|^_,, 

where, according to the notation (II. lip . stands for 

5„n|_, = -e, • (1 + QH)^^'^^.,. 

Before stating the main result of this section let us introduce a nota- 
tion: 

Notation 2.1. We generically write 

M[a] := C{e^, ||V^'^^7||iooj^*o+i) , (2-9) 

where, as usual, C(-) is a nondecreasing function of its arguments. 

Proposition 2.4. Let to > 1, s > and C,b e H^°+'^ n i7*+i(M2) be 
such that (|2.4p is satisfied, and let a be given by (|2.5p . 
Then for ah g G C{[-1,0]; H''{R'^f), there exists a unique variational 
solution u £ H^{S) to the EVP and 

P^V^'^TxIb <MH(||.l^g||2 + (||V'^'^a||L^^,.||.l*°g||2>,>,^+J, 

where M[a] is defined in (j2.9p . 



Proof. The existence of the solution can be obtained with very clas- 
sical tools and we therefore omit it. We thus focus our attention on 
the proof of the estimate. 

Let x(") be a smooth, compactly supported function such that xiO — 
1 in a neighborhood of ^ = 0, and define := A*x{^D). Using Af^^u 
as test function in the variational formulation of (|2.8p . one gets 

[ (1 + Q[(t])V^'^u- V^'^yl^^n = [ g-V^'^4^M, 
Js Js 

so that using the fact that Af^ is L^-self-adjoint, one gets, with Vh = 
Alu, 



[ AUl + Q[a])V^'''u-V>^'^VH= [ Aig-V^- 
Js Js 



and thus 



Js Js 
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Thanks to the coercitivity property of Proposition 12.31 one gets 

k[a]-^\\AiV'^''^uh < \\[AlQ[a]]V^M\2 + \\Aig\\2; (2.10) 

since is of order s — 1, the above estimates allows one to 

conclude, after letting h go to zero, that A^V^^'^u G Lp'{S); more 
precisely, thanks to Corollary 12.11 oris deduces 

and thus 

WA'V^^^uh < C{k[a], ||gM|L^^.o+i) (2-11) 

X {\\A^g\\2 + \\V^'M\2 + {\\Q[a]U^H4^'''V^'M\2)s>to+i)- 

One also gets ||V^''^n||2 < A:[(T]||g||2 from (I2.10p after remarking that 
the commutator in the r.h.s. vanishes when s = h = 0; taking s = to 
in (EH]) then gives V^'^^lb < C{k[a],\\Q[a]\\^^^,^+,)\\A'^'g\\2, 
so that the r.h.s. of (j2.1ip is bounded from above by 

C(M^],||QMIlLoo//*0+i)(P^g||2 + (||QH||L^/^.||yl*«g||2),>,^J. 

The proposition follows therefore from Proposition 12. 3i □ 

Before stating a corollary to Proposition 12. 4^ let us introduce a few 
notations: 

Notation 2.2. i. For all u G H^/'^(M?), we define as the solution 
to the BVP 

fv/^'^-(l + Q[a])V'^'V = 

ii. For all u G ©'(R^), one defines as 

yze[-l,0],u^;z)=x{VJIz\D^\)u, 

where x is a smooth, compactly supported function such that x(0) = 
1. 

The following corollary gives some control on the extension mapping 
u v!' . 

Corollary 2.2. Let to > 1 and s > 0. Let also C,b £ i7*o+2 n 
^s+i(^]^2^ be such that ([23]) is satisfied, and a be given by (f23D . 
Then for all u G H^'^^^'^{M?), there exists a unique solution G 
H^{S) and 

with *P as defined in 
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Proof. Looking for under the form v!' = v + u\ with given by 
Notation 12.21 one must solve 



V^'^ • (1 + Q[a])\/''^^v = -V^'T • (1 + gi^DV^'^nt 



(2.13) 



, v\,=o = 0, dnv\,=.i = e, • (1 + QH)V^'VU=_i. 
Applying Proposition 12.41 (with g = — (1 + Q[a])V^''^u''), one gets 

and since = + v, the corollary follows from Proposition 12. 2i □ 

Remark 2.2. From the variational formulation of (j2.13p . one gets 
easily 

11(1 + QM)l/2yM,7^||^ <\\il+ gH)l/2yM,7^t||2. 

3. The Dirichlet-Neumann operator 

As seen in the introduction, we define the Dirichlet-Neumann opera- 
tor Gfj,,^[eC,[3b]- as 



where <P solves (j2.3p . Using Notation 12.21 one can give an alternate 
definition of (see Proposition 3.4 of [29j), namely, 

g^,^[eC,/3&]^ = 5„/U=o ( = • P[a]Vx,./U=o)- 
More, precisely one has: 

Proposition 3.1. Let to > 1, s > and C, & e H^°^^ n i?^+^(M2) be 
such that (j2.4|) is satisfied, and let a be given by (j2.5|) . 
Then one can define the mapping Q^^-yleC, (3b]- (or simply ^[eC]" when 
no confusion is possible) as 

G,^,yH,Pb]{=g[eC]) : „ ' fl „b| • 

U OnU \z=0 

Proof. The extension is well-defined owing to Corollary 12.21 More- 
over, we can use the definition of P[(t] and Qlc] to see that 

ez • PHVx.^n^ = ez • (1 + Q[a])V''^^u\ 

We will now show that it makes sense to take the trace of the above 
expression at z = 0. This is trivially true for Q[a], so that we are left 
with . After a brief look at the proof of Corollarv 12.21 ^^'i using 
the same notations, one gets = v + u"^ . Since one obviously has 



16 



B. Alvarez-Samaniego, D. Lannes 



V^''^v) G C{[-l,Q]]H''-^/'^{m?f), the trace V^''^v)\z=Q makes sense. 
In order to prove that V^'^v\z=q is also defined, remark that f , which 
solves (l2l3D satisfies V^'^w G L2((-1, 0); i?'(M2)3) and, using the 
equation, dzV''^v G L^{{-1,0); H'^-^ (M.'^f). By the trace theorem, 
these two properties show that V^-'^'f 1^=0 G H^~^/'^(M?)^ . □ 

3.1. Fundamental properties 

We begin this section with two basic properties of the Dirichlet- 
Neumann operator which play a key role in the energy estimates. 

Proposition 3.2. Let to > 1 and C, 6 G F*o+2(]^2) ^^^^i that (pli]) 
is satisfied. Then 

i. The Dirichlet-Neumann operator is self-adjoint: 

Vn, V G fV2(m2), (u, g[eC]v) = {v, g[eC]u). 

ii. One has 

yu,v G H^/\R^), \{u,g[eC]v)\ < {u,g[eC]u)^/\v,g[eC]v)^/^ . 
Proof. Using Notation 12.21 one gets by Green's identity that 

{u,g[eC]v) = /(I + Q[a])V^''"'u^ ■ V^'V (3.1) 
Js 

= j(l + QM)i/2v'^'^u^ . (1 + QH) V2 V, (3.2) 

where (1 + Q[c])^'^^ stands for the square root of the positive definite 
matrix (1 + Q[o"]) (note that the symmetry in u and v of the above 
expression proves the -very classical- first point of the proposition). 
It follows therefore from Cauchy- Schwartz inequality that 

{u,gH]v) < 11(1 + QH)i/2v'^'V||2 11(1 + QH)i/2v^'V||2, 

which yields the second point of the proposition, since one has 

iu,g[£C]u) = ||(l + QM)i/2yW||2 (3.3) 

(just take u = v in (j3.2p ). □ 

The next proposition is related to the variational formula of Hadamard 
and gives a uniform control of the operator norm of the DN operator 
and its derivatives (recall that we use the convention (jl.lOp and that 
dig[e-]u ■ h = g[e(]u when j = 0). 
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Proposition 3.3. Let to > I, s > and C,b £ H^o+'^ n H''+^{R'^) be 

such that (|2.4p is satisfied, and let a be given by (|2.5p . 

For aU u e H'+^/'^{R'^), j E N and h G p H'+'^{M.y, one has 

1 / 
I— 4gHn- h|^,_i/2 < (-yMH(|*Pn|^, n l^-tlj/^o+i 

j 

+((1 + IIV^'VIIl^hOI^P^Ih^o n \^k\mo+^),^to 

k=l 

j 

IT I IT I 



5>tO 

fe=l /^fc 

with M[a] as in ([23]) while *p is defined in (fL8]) . 

Remark 3.1. When j = 0, the proposition gives a much more precise 
estimate on |^[eC]^l_ff''-i/2 than Theorem 3.6 of [2^, but requires ( G 
while C G Z^'+^/^lM^) is enough, as shown in [29] through 
the use of regularizing diffeomorphisms. This lack of optimality in 
the ("-dependence is the price to pay to obtain uniform estimates in 
terms of a <S^(-) rather than Sobolev-type norm. 

Remark 3.2. The r.h.s. of the estimate given in the proposition 
(when j = 0) is itself bounded from above by 

M[a]{\u\Hs+i + (l|V^'V||L-//H^^|H*o+i>s>to)- 

Proof. First remark that one has A'^^^/'^v'^\z=o = A'^~^/'^v (with f''' as 
in Notation 12. 2p . so that one gets by Green's identity. 



iA'-'/^g[eC]u,v) = {g[eC]u,A'-'/^v) 

■u'. (3.4) 

A Cauchy- Schwartz inequality then yields, 

iA'-^/^g[eC]u,v) < ||yl^(l + QH)V^'V||2M-^/2v'^'V||2, (3.5) 

and since it follows from the product estimate (|2.ip that ||vl''(l + 
Q[cr])V^^''^u^\\2 is bounded from above by 



(1 + ||QH||^^^*o)P^vW||2 + (||QH||ioo^,.p*ovMV||2) 



S>to' 
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one can deduce from Propositions 12.21 and 12.31 that (recall that v = 

and the proposition thus follows directly from Corollary 12.21 and a 
duality argument in the case j = 0- 

In the case j ^ 0, after differentiating (j3.4p . and using the same 
notation as in Proposition 12.31 one gets 

{A'-^/^dig[e-]u ■h,v)= I A\Q^^\a] ■ h)V'^'^n^ • A-^/^V^'^^v^ 

J s 

+ j A'B{hk,\y,-k)-A-^/^V^'^v\ (3.6) 

k=l hfe,hj_fc ^ 

where the second summation is over all the /c-uplets and (j — k)- 
uplets h.j_k such that {hk,h.j_k) is a permutation of h, and where 
B{h.k,h.j-k) is given by 

B(h,,h,_fe) = (Q(^'-^)M ■ h,_fe)V^'^(n^''= • hfc) 

{v}''^ ■ h^. standing for the k-ih. order derivative oi Q ^ v!' a.t C, and in 
the direction h^). 

Proceeding as for the case j = and using the estimates on 1 1(5 ['^] " 
h||Loo_f/s provided by Proposition 12.31 oris arrives at the desired esti- 
mate for the first term of the r.h.s. of (j3.6p . For the other terms, one 
has to remark first that v}''^ ■ solves a bvp like (j2.8p with 

g = -X] Y B{hk,i,hk^k~l), 

1 = hfc,;,hfc,fe_; 

where the second summation is taken over all the / and k — /-uplets 
such that (hfc h^^^.;) is a permutation of h^. A control of \\V^'^v!''^ ■ 

hfc||L°°H= in terms of \\Q^''~'-'>[a]-hk,k-i\\L^H'' and \\V^'"'u^''' ■hk,i\\L°°H'' 
is therefore provided by Proposition 12. 4[ It is then easy to conclude 
by a simple induction. □ 

Remark 3.3. Instead of (|3.5p . one can easily get 

{A'-^/^g[eC]u,v) < \\A'+'^/'^{l + Q[a])V>''"'u^\\2\\A-^V^''"'v^\\2, 
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and since ^V'^''^u^||2 ^ y^|w|2 one also has the estimate (with 
C = for the sake of simphcity) : 

1 1 |DT| 

showing that -^[Ol- can be uniformly controlled when /i goes to zero. 

The proposition below show that controls in terms of |^m|2 or 
(n, -^Q[£C\'^Y^'^ ^-re equivalent. This result can be seen as a version 
of the Garding inequality for the DN operator. 

Proposition 3.4. Let to > 1 and C,^ G i7*o+2(]R2) be such that 
()2.4p is satisfied, and let a be given by (|2.5p . A;[cr] be as defined in 
Proposition O and ^ be given by For all u G H^I'^{'E?), one 

has 

{u, —G[eC,]u) < M[a]\^u\l and k[a]-^\^iu\l < (n, —g[eC]u). 



Proof. The first estimate of the proposition follows directly from 
and Corollarv 12.21 

The second estimate is more delicate. Let (/9 be a smooth function, 
with compact support in (—1,0] and such that Lp(fi) = 1; define also 
u(X, z) = Lp{z)v!' (with defined as in Notation 12. 2p . Since v\^^_^ = 
0, one can get, after taking the Fourier transform with respect to the 
horizontal variables, 

|£7|2 /-O |£7|2 

Remarking that 

l^^l < Iv^lool^^l and \dzv\ < \dz(p\oo\u^\ + \<f\oo\dzU^\, 
one gets 

|t7|2 rO |t7|2 _ 



rO ,,|C7|4 ^ /-O 1 

|2 / /^l? I L.b/f .m2,„ , I, i2 / ^\a..\,/t „m2. 



+l'^l°o / (i + ^|^7|)2 l^^(^'^)r^^ + l'^ISo j_^-\dzu^{^,z)fdz, 
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where Young's inequahty has been used to obtain the last hne. Re- 



marking now that ^ /^j^iLuij < IC^P, one has 



so that, integrating with respect to ^, one gets 

^ < -||V^'V||2 



\D^\ 



2 /X 



(l + V^|i)7|)l/2 

Owing to Proposition 12.31 and ()3.ip (with v = u), one has 

||V^'V||2<A:M(n,gK]n), 
and the proposition fohows. □ 

3.2. Commutator estimates 

In the following proposition, we show how to control in the energy 
estimates, the terms involving commutators between the Dirichlet- 
Neumann operator and spatial or time derivatives and in terms of 
rather than Sobolev-type norms. 

Proposition 3.5. Let to > 1, s > and C,b e H'^"^^ n H^^^iR"^) be 
such that (12. 4p is satisfied, and let a be given by (12. 5p . Then, for all 

V G 7/^+1/2 (M2), 

\[A',^g[eC]]v\^<M[a](^\\V^''^ah.cHt,+i\^v\Hs 

where M[a] is as in (j2.9p while %i is defined in (|1.8p . 
Proof. First remark that for all u G ©(M^)^ 

(u, [g[eC],A']v) = {u,g[eCWv) - {A'u,g[eC]v) . 

Since u^|2=o = u and A'^u^\z=o = A^u (we use here Notation 12. 2p . it 
follows from Green's identity that 

(n, [g[£C],A']v) = V^'^nt • (1 + Q[cj])V^'^'< {A'vf 

- / (1 + Q[a])V>''^v'' ■ V^'^^iA'v) 
Js 

= j V^'V • {{I + Q[a])V^'^^ {{A'vi' - A'v^) - Q^IV^V). 

(3.7) 

Let us now prove the following lemma: 
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Lemma 3.1. For all / G ^^(M^) and g G H^{Sf, one has 

J^y^''P-g<^^V^\f\2\\Ag\\,. 
Proof. By definition of p , one has 

Replacing V'^''''/^ in the integral to control by this expression, and 
using the self-adjointness of A, one gets easily from Proposition 12.21 
that 

Js 

recalling that u = a-nd 7 < 1, one can check that |*pyl~"'^/|2 ^ 

I/I2, uniformly with respect to /j, and 7, and the lemma follows. □ 



It is then a simple consequence of the lemma, (|3.7p and Proposition 
[211 that 



{u, [g[eC],A^]v) < ^V^\u\2[\\A[A^Q[a]M^^v'h 

+ (1 + ||V^'V|L^^.o+i)PV'^'"((yl^^)^ - rv')h), (3.8) 

which motivates the following lemma: 
Lemma 3.2. One has 

\\AVf''"'{{A'vf - A'v^)\\2 < M[a]\\A[A' ,Q[a]]Vf''"'v^\\2. 

Proof. Just remark that w := (A'^v)^ — A^v^ solves 

V^'^ • (1 + Q[a])V''^^w = V^'^ • g, 
w\^=o = 0, 9„w;|^=_i = -e^ • gU=-i, 

with g = [A^, Q[a]]V^''^v^ , and use Proposition [231 □ 
With the help of the lemma, one deduces from (13. 8p that 

{u,[g[eC],A']v) < ^V^M[a]\\A[A',Q[a]]V>'^^v^\\^\u\2, 
and thus, owing to Corollarv 12.11 and Proposition 12.31 

{u,[g[eC],A']v) < ^V^M[a]\u\2 

and the result follows therefore from Corollarv 12.21 and a duality ar- 
gument. □ 
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The next proposition gives control of the commutator between the 
Dirichlet-Neumann operator and a time derivative. 

Proposition 3.6. Let to > 1, T > and C,b e C^{[0,T]; H^o+^{'R'^)) 
be such that ()2.4p is satisfied (uniformly with respect to t), and let a 
be given by ([23]). Then, for all u G C^{[0, T]; H'^/'^{M.'^)) and t G [0, T], 

\{[dt,^g[eCMt),u{t))\< M[ait)] \\V^''Wta\\^,T\V<t)\l 

where M[a{t)] is as in ([^ while ^ is defined in piH]) . 
Proof. First remark that 

(^i, [at,gK]]t^) = at(n,gK]n) - 2{u,g[eC]dtu), 
so that using Green's identity, one gets 



-2 y (l + QH)V'^'^(9tn)^- V^'V 

-2^(1 + Q[a])V'''^{idtuf - dtu^) ■ V^'^Tx^ 

It follows directly that 

{u,[dt,g[eC]]u) <\\dtQ[a]\U\V^'"'v}'\\l 

+ (1 + WQMUWV^'^iidtuf - atn^)||2||V^'V||2. 

Proceeding exactly as in the proof of Lemma 13.21 one gets 

WV^^'^idtuf -dtu^)\\2<\\dtQ[a]\U\Vf''"'u^\\2, 
and the result follows therefore from Corollary 12.21 and Proposition 

E31 □ 



3. 3. Other properties 

Propositions 13.2] and [3^ allow one to control {u,Q[eC]v) in general. 
However, it is sometimes necessary to have more precise estimates, 
when u and v have some special structure that can be exploited. 
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Proposition 3.7. Let to > I, s > and C,b £ H^o+'^ n H''+^{R'^) be 
such that (|2.4p is satisfied, and let a be given by (|2.5p . 

i. For all v G H'+^ n F*o+2(m2)2 ^ ^ H'+'^/'^{M.'^), one has 

{[A',y] ■ V^u, ^e?[eC]([yl^ v] • V^n))'/' < M[a] 

ii. For all v e iJ'o+VM^)^ and u G i7i/2m2N j^^^g 



((v • V^u), ^G[eC]u} < M[a] \x\w^,^ l^^ug 

Proof. In order to prove the first point of the proposition, define 
U = [^*, v] • V^u and recall that we saw in p.3p that 



{U,g[e(:]U) = ||(l+g[c7])l/2yM,7[/b|| 



Using Notation \2.'2\ we define = [yl^v] • V''' u"^] as in Remark 
(with instead of t/t), one deduces 

{U,G[eC]U) < \\{1 + Q[a])^/^V^''^U%. 

Since 7 < 1, one has \\V^U% < \\AU% and one gets with Proposi- 
tion Ell 

llV^C/^lb < |v|^.o+iM^vV||2 + (|vU.+i||.l*«vV||2>,>,^; 
similarly, since dzU'^ = [^'', v] • V^S^n^, one gets 

\m% < |v|^*o+iP^5,nt||2 + (|v|hH|v1*°+^9,J||2>,>,„+i, 
so that, finally, 

l|V^'^C/^||2<|v|^*o+2P^V^'V||2 + (|vU.+i||yl*«+iV^'V||2),>,^^,. 
It follows therefore from Proposition 12.21 that 



l|V^'^C/^||2<^|v|^*o+2l 



\D'^\ 

+ (V^|v|hs+i| ^ ^|D7|)l/2 "lg'o+^)^>to+l' (3-9) 

and the result follows. 

To establish the second point of the proposition, first remark that, 
owing to Green's identity, 

((v • V^n), g[eC]u) = / (1 + Qi^DV^'Tn^ • V^'^(v • V^n^), (3.10) 
Js 
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SO that, 

{{^■V"'u),g[eC]u) = J (1 + Q[(7])V'^'V- [V'^'^,v-V^]n^ 
+ J V^^'^'u^ ■ [Q[a], (v • V^)]V^'^u^ 

+ / V^'V- (v.V^)(l + QH)V^'^u^ (3.11) 
Js 

Integrating by parts, one finds 

J V^'V-(v-V^)(l + QH)V^'V 

= - / ((div^v) + V • V^) V^'^n^ • (1 + Q[a])V^''^u^ (3.12) 
Js 

= - J (div^v)V^'^M^ • (1 + Q[a])V^''^v!' 

- y [v • v^, ■ (1 + QM)v^'^u^ 

-J ^^'^(v- VV) • (l + Q[a])V'^'^u^ (3.13) 

From (IXTTIl . (IXT2I1 and (f3l3]) . one gets therefore 

((v-V^n),g[eC]^x) = y (l + Q[(j])V'^'^u^-[V^'^,v-V^]u^ 
+ i /" v'^'Tu^ • [QH, (v • V'^)]V'^'^n^ 

-If (div^v)V^'^'u^ • (1 + QM)V^'^'«^ 
2 Js 

Remarking that [V^'^, v • V^] = + V^VsT^^^j,^ ^ 

deduces easily that 

((v • V^n),g[eC]u) < |v|h^i,oo(1 + IIQHII^i.oo)!!^^'^!!^ 

and the results follows from Corollary 12. 2i □ 

We finally state the following theorem, which gives an explicit 
formula for the shape derivative of the Dirichlet-Neumann operator. 
This theorem is a particular case of Theorem 3.20 of 
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Theorem 3.1. Let to > 1, s > to and C,b e H^+^/'^(E?) be such that 
(123D is satisfied. For ah 'ip e H''+''^/'^{R'^), the mapping 

is well defined and differentiable in a neig hborhood of ( in (1^2 ^ ^ 

and 

V/i G /7^+3/2(m2)^ d^g[e-]i^ ■ h = -eg[eC]{hZ) - e^V^ • (/iv), 
with Z := 2^[eC]V' and v := V^ip — sZ\/^(, and where 

Remark 3.4. We take this opportunity to correct a harmless mis- 
print in the statement of Theorem 3.20 of [29j. It should read 



daG{;b)f-h = -G{a,b){hZ)- 




and Pa should be replaced by P on the right hand side of the equation 
in the statement of Lemma 3.24. 



3.4- Asymptotic expansions 

This subsection is devoted to the asymptotic expansion of the DN 
operator G[e(]'4'{= ^^^^[eC; iii terms of one or several of the 
parameters e, fj,, 7 and p. We consider two cases which cover all the 
physical regimes described in the introduction. 

3.4-1- Expansions in shallow-water <^ 1) In shallow water, that 
is when /i <C 1, the Laplace equation (|2.3|) -or its straightened ver- 
sion (|2.6P " reduces at first order to the ODE = 0. This fact 
can be exploited to find an approximate solution ^app of the Laplace 
equation by a standard BKW expansion. This method has been used 
in the long- waves regime in [5l[32l[9] (see also [39]) where the corre- 
sponding expansions of the DN operator can be found. We prove here 
that it can be used uniformly with respect to e and /?, which allows 
one to consider at once the shallow- water /Green-Naghdi and Serre 
scalings. The difference between both regimes is that e = /? = 1 in 
the former (large amplitude for the surface and bottom variations), 
while e = P = ^/JI in the latter (medium amplitude variations for the 
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surface and bottom variations). 

Let us first define the first order linear operator T[/i, b] as 

T[h, b]V:=-^V{h^V ■ V) + \[Vih^Vb ■ V)-h^VbV ■ V] +hVbVb ■ V. 

(3.14) 

Proposition 3.8 (Shallow- water and Serre scalings). Let 7 = 

1, s > to > 1, G -H'*+"/^(IR^), b G i?*+ii/2(]K2) ^ g 
^s+9/2|-]^2) assume that ([23]) is satisfied. 
With /i := 1 + — one then has 

with Cj = C{-^,\C\H^+5/2+2j,\b\jjs+r/2+2j,\V'^tp\Hs+7/2+2j) {j = 0,1), 

and uniformly with respect to e,/? G [0, 1]. 

Proof. We look for an approximate solution (papp to the exact solution 
(p of the potential equation (j2.6p under the form 

4>app{X, Z) = lp{X) + flCpiiX, z). 

Plugging this ansatz into (j2.6p . and expanding the result into powers 
of fi, one can cancel the leading term by a good choice of namely, 

z^ 

MX,z) = -h{h{— + z)A^ - zf5Vb -Vi)). 

One can then check that 

^x,z ■ P[cr]'^x,z4'app = fJ-'^R^,, in 5, 

with (R^jr/j^) satisfying, uniformly with respect to /x G (0, 1), 

r^\jjs + l/2 < C{\C\f^s+5/2,\b\jjs+7/2,\V^'lp\Hs+7/2). 

(3.15) 

Smce g[eC]ip - dncpappi^^o = 5„(0 - (/'app)|^=o' the truncation error 
can be estimated using the trace theorem and an elliptic estimate on 
the BVP solved by (p 

~ (papp'i this is exactly what is done in Theorem 
1.6 of [9j for instance, which gives here: 

\G[eC]^ - dn(t>app uJh^ < ^l''Cs{\W^^'^R^\\2 + |'"/i|//s+l/2), 

with Cs = C{\C,\fjs+i,/2, \b\fjs+i,/2). Together with (j3.15p . this gives the 
result. 

In order to prove the second estimate of the proposition, one must 
look for a higher order approximate solution of ()2.6p . namely (papp = 
Tp + ficpi + /i^</>2. The computations can be performed by any software 
of symbolic calculus and the estimates are exactly the same as above; 
we thus omit this technical step □ 
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3.4-2. The case of small amplitude waves (e 1) Expansions of 
the Dirichlet-Neumann operator for small amplitude waves has been 
developed in |17pi6j. This method is very efficient to compute the 
formal expansion, but instead of adapting it in the present case to 
give uniform estimates on the truncation error, we rather propose a 
very simple method based on Theorem 13.11 

Proposition 3.9. Let s > to > 1, W S H^+'^/'^{^^) and C G 

^s+3/2(^]^2^) g^p^ ^^YysX (fTIj) is satisfied for some /iq > 0. Then 
one has 

U tlQ 

Proof. A second order Taylor expansion of Q[eQip gives 

gK]^ = a[0]V + doa[£# • C + /V - ^)dlS^-]^ ■ (C, C)dz. 

Jo 

Using Theorem 13. H one computes 

doa[e# • C = -eg[0]iCiGm)) - efiV^ ■ (CV^V), 

while for all z G [—1,0], Proposition 13.31 controls dl^^Qle-jilj ■ (CiC) iii 
-ff* by the r.h.s. of the estimate given in the statement. □ 



(4.1) 



4. Linear analysis 

The water-waves equations (|1.4p can be written in condensed form 
as 

dtU + CU + -A[U] = 0, 

V 

with JJ = {C,i'V, AU] = iAi[U],A2[U])'^ and where 

-1 
and 

Ai[u] = -j-^{g[eC]i^-gm, 

1 (^gKl^ + sy^v^C-v^V^)^ (4.2) 

MU] = -\V^i^\ 2(l + eV|V7CP) • 

By definition, the linearized operator around some reference 

state U_ = (C, ip)'^ is given by 

= dt + C + ^duA; 
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assuming that U is such that the assumptions of Theorem 13.11 are 
satisfied, one computes that ^{(,^) is equal to 



dt + 



-2 



(4.3) 

where v and Z are as in the statement of Theorem 13.11 

This section is devoted to the proof of energy estimates for the 
associated initial value problem, 



Defining 



a = 1 + ^b, and b = ev • V^Z + udtZ, (4.5) 



we first introduce the notion of admissible reference state: 

Definition 4.1. Let to > I, T > and b e H^o+^{R'^). We say that 
U_ = (C, ip) is admissible on [0, ^] if 

— The surface and bottom parameterizations C and b satisfy ()2.4p 
for some /iq > 0, uniformly on [0, 

— There exists cq > such that a > cq, uniformly on [0, 

We also need to define some functional spaces and notations linked 
to the energy ()1.6p mentioned in the introduction. 

Definition 4.2. For all s G M and T > 0, 

i. We denote by the vector space H'(R'^) X if^+i/2(M2) endowed 
with the norm 

yu = {c,4^f ex', \u\xs -IClH^ + ^mHs + miH^, 

while stands for C([0, ^];X^) endowed with its canonical norm. 

ii. We define the space X^ as 

:= {U = (C,'0)^,C e H'{R'^),Vil^ G H'-^/\R^f}, 

and endow it with the semi- norm := ICIh" + l^^l/f^- 

iii. We define the semi-normed space {Y^, \ ■ |y^) as 

YS:=^CH[0,'^];X^-h and |t/|y^ =^ sup l^.^^f/l s,. 
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iv. For all (G, G X|, x X^ we define 

T\t,U\G):=\U^\xs + - f sup \G{t")\xsdt' . 

V Jo 0<t"<t' 

We can now state the energy estimate associated to (14. 4p . and whose 
proof is given in the next two subsections. 

Proposition 4.1. Let s > to > 1, ^ > 0, 6 G H''-^^/'^{R'^), and 
U_ = (C, "0) G Y^~^^^'^ be admissible on [0, ^] for some /iq > and 
CO > oTh^t also (G, [/°) G ^ ^s+2_ 

There exists a unique solution U G to (I4.4jl : moreover, for all 
< t < ^, one has 

where C = G(r, ^, f , f , |6|h=+9/2 , iC/L^o+s/^) . 



^.i. Energy estimates for the trigonalized linearized operator 

As shown in [29], the operator is non-strictly hyperbolic, in 

the sense that its principal symbol has a double purely imaginary 
eigenvalue, with a nontrivial Jordan block. It was shown in Prop. 
4.2 of [29] that a simple change of basis can be used to put the 
principal symbol of ^((^ip) under a canonical trigonal form. This result 
is generalized to the present case. More precisely, with as defined 
in (j4.5l) and defining the operator = dt + with 

one reduces the study of (j4.4p to the study of the initial value problem 

as shown in the following proposition (whose proof relies on simple 
computations and is omitted). 

Proposition 4.2. The following two assertions are equivalent: 

- The pair U = (C, -0)^ solves 

-The pair V = {C,^ - eZCf solves gZ]), with H = (Gi,G2 - 
eZGif and = (C°, V'^ - e^|,=oC°)^- 
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In view of this proposition, it is a key step to understand 
and the rest of this subsection is thus devoted to the proof of energy 
estimates for this initial value problem. 

First remark that a symmetrizer for is given by 

S=(nsl^lr.i.,^.)^ (4-8) 



so that (provided that a is nonnegative), a natural energy for the IVP 
(j4.7p is given by 

E%V)'^ = {A'V,SA'V) 

= \V^^'V,\l + ^1^21^. + iA'V2, ^g[eC]A'V2). (4.9) 



Remark 4.1. The introduction of the term e'^/v'^ in (j4.8p -and thus 
of e^/j/^|V2||/s in (|4.9|) " is not necessary to the energy estimate below. 
But this constant term plays a crucial role in the iterative scheme used 
to solve the nonlinear problem because it controls the low frequencies. 
It also turns out that the order 0(e^/z^^) of this constant term is the 
only one which allows uniform estimates. 

We can now give the energy estimate associated to (|4.7p : in the 
statement below, we use the notation 

r{t,V\H):=E'{V^) + - f sup E'{H{t"))dt', 

^ Jo 0<t"<t' 

while s\/ to := max{s, to} smd C_ is as defined in Proposition 14. 1[ 

Proposition 4.3. Let s > 0, to > I, T > 0, b e H^'^^o+^/'^{R'^), and 
U_ = (C, V') £ 5/^^*0+9/2 admissible on [0, — ] for some /iq > and 
Co > OT ~ 

Then, for all {H,V^) G x X^, there exists a unique solution 
V £ X^to (132D and for ah < t < ^, 

E'{v{t)) < c{r{t,v',H) + (|f/|^^+./.i*«+i(t,yo,i/)>^^^^^J. 

Proof. Throughout this proof, C_q denotes a nondecreasing function 
'^f Co' ^' "^^[— ]' IxIh*o+2) |J^|j:/to+2, and \dtb\oo which may vary from 
one line to another, and a is given by 12. 5\) with C = C- 
Existence of a solution to the IVP (|4.7p is achieved by classical means, 
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and we thus focus our attention on the proof of the energy estimate. 
For any given k G M, we compute 

-2{A'M(^^^^)V, SA'V) + {A'V, [dt, S]A'V). (4.10) 

We now turn to bound from above the different components of the 
r.h.s. of (glOj) . 

• Estimate of {A^H, SA'^V). We can rewrite this term as 

iV^A'Hu ^A'Vi) + {-A'H2, -A'V2) + {A'H2, —g[eC]A'V2), 

so that Cauchy-Schwartz inequahty and Proposition 13.21 yield 

{A'H, SA'V) < E'{H)E'{V). (4.11) 

• Estimate of {A^M(^(^^^-^V,SAW). One computes 

(^U%,^)y,S^^y) = (^^(Jdiv^(vyi) - l-g[eQV2),aA'Vi) 

+ {A^{aVi + -V • V^Fs), (4 + —Q[eQ)A^V2). 

so that one can write 

(^^A'%^)y, SA'V) = ^1 + yl2 + ^3 + ^4 + ^5, 

with 

Ax = -(A'divJvyi),aA'Vi), 
1/ ' 

A2 = -[A\^-V^V2),t^A'V2). 

^3 = -(yl^(v. V^Va),— a[eC]^'V^2), 

^4 = {A\^Yx),—Q{eC,\A'V2) - (oyl^Fi, yl^(— aKl^s)), 

A^ = [A\aV{),t^A'V2). 
We now turn to prove the following estimates: 

A, < ^C,E%V)[{1 + pV'^'^a||^^^.o+.)i?^(y) 

+{{\y\hs+^ + + pV'^'^a||io.^.+i)ii;*»+H^)>,>i„+i),(4.12) 

for j = 1, . . . ,5. 



32 



B. Alvarez-Samaniego, D. Lannes 



— Control of Ai and A2- Integrating by parts, one obtains 

-Ai = ([^^div^(v•)]14,ayl^^l) - l{A'Vi,{v-V^a)A'Vi) 
+\{aA'Vi,{diw^Y)A'Vi) 

and 

2 2 

-^A2 = ([yl^v]V^y2,^/l^■^^2) - ^{A'V2,{div^^)A'V2). 

Recalling that a = 1 + one can then deduce easily (with the 
help of Proposition 12.11 and Corollary 12.11 to control the commu- 
tators in the above expressions) that (j4.12p holds for j = 1, 2. 

— Control of A^,. First write ^3 = ^31 + ^32 with 

^31 = -{ — QHWV2, V] • V^V2) 

Az2 = -(v • V^yl^Fa, — a[eC]^'^2) ■ 

I' — 

Thanks to Proposition 13.21 one gets 

A31 < v] • v^y2, v] • v^'V2Y'^E\v), 

and Propositions I3.7l i and 13.41 can then be used to show that A31 
is bounded from above by the r.h.s. of (j4.12p . This is also the case 
of ^32, as a direct consequence of Propositions I3.7l ii and 13. 4i It 
follows that (fiT2|) holds for j = 3. 

— Control of ^4. One computes, remarking that [^^,0] = ^[yl'^, b], 

A4 = {aA''Vi,[—QH],A'^]V2) + -{[A\h]Vi.—QHWV2) 
:=^4i +^42. 

Using successively Cauchy-Schwartz inequality, Proposition 13.51 
and Proposition 13.41 one obtains directly that ^441 is bounded 
from above by the r.h.s. of (|4.12|) . In order to control ^42, first 
remark that using Propositions 13.21 and 13.41 one gets 

^^42 < {[A\b]Vi,^Q[eQ[A^MyiY'^E^{V) 
<M[a]\^[A\b]Vi\^E\V). 
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Recalling that v = jif^ one can check that for all ^ G M^, 

(i^^|g7|)i/2 ^ (0) uniformly with respect to n and 7, so that 



-l/2|g7 

one deduces 

^Ai2 < M[a]\[A',b]Vi\^,E'{V). 

Remarking that owing to Proposition 12. H one has 



and A42 is thus bounded from above by the r.h.s. of (j4.12p . This 
shows that (j4.12p also holds for j = 4. 
Control of A^. First remark that 



A5 = {A'Vi, ^A'V2) + -{A'ibVi), ^A'V2), 

so that Cauchy-Schwartz inequality and the tame product esti- 
mate dZj]) yield 



^5 <-((! + \-b\H^,)\yiW + (i-^iHH^ib*o>,>.,)-i^2k» 

< "-^^{(^ + \"-b\H^o.^)E%V) + {\'-b\H^E''^^\V))^^^^^^)E%V), 

and (|^TT^ thus holds for j = 5. 
From (|4.12|) . we obtain directly 

{A^M^^_^^)V,SA^V) < ^CoE'^{V){{l + ^||V^'V||io.^to+2)i?^(y) 

(4.13) 

• Estimate of {A'V, [dt, S]A'V). One has 

(A'V, [dt, S]A'V) = ^{A'Vi, dtbA'Vi) + {A'V2, [du ^g[eCWV2), 

so that, using Proposition 13.61 to control the second component of the 
r.h.s., one gets easily 

(A'V, [dt, S]A'V) < JCo(l + "^WV^'^dtalUE^Vf. (4.14) 
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According to (fiJOl) . (filT]) . KTSh and we have 

(4.15) 

with := (|v|j^s+i + |||V^''>'a||j;^ooj|^s+i + l^l^^+i)' Provided that k 
is large enough, how large depending only on 



sup 



Co(t)(l + -||V'^'Mi)llLooj,*o+2 + -||V'^'^9ia(t)||, 
It follows from (|4.15p that, 

V Jo 

+ {-Coi sup Ds) [ ev(*-*')^to+i(y(t/))^i'\ 

lO.vT/e] Jo 



's>to+l' 

(4.16) 

using (|4.16p with s = to + 1 gives 

£;*o+i(y(t)) <eV*^*o+i(yO)^£^g^t sup E'^'+\H{t')), 

V 0<t'<t 

and plugging this expression back into (I4.16P gives therefore 
E%V{t))<C,{l\t,V',H) + {{ sup Z),)/*°+i(t,F0,i7)> ), 

te[0,uT/e] 

where C_i is a nondecreasing function of T, ^ , ^ , ^ and of the supre- 

mum on the time interval [0, ^] of f || V'^''^a||^cx)j|^to+2, ^\\V^''^dta\\oo, 
1x1^*0+2, |b|j|^to+2 and |5tb|Lcx). The proposition follows therefore from 
the following lemma: 

Lemma 4.1. With C_ and | • as defined in the statement of Propo- 
sition HTT] and Definition 14.21 one has, 

Vs > to + 1, sup Ds{t) < C\U\ s+7/2 and < C. 

telO,uT/e] T 

Proof. Remark first that, as a consequence of Proposition 13.31 
has for all r >tQ + 1, 



(4.17) 
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since < ^ (i+vT^I^"^) ^ ' uniformly with respect to 7 G (0, 1] and 
fj, > 0, one also has 

\V^t\Hr < + < \U\xr+^/2- (4-18) 

It follows from the explicit expression of Z given in Theorem 13.11 that 
eZ_ is a smooth function of e^/Jl < e/v, VC, V^ip and -^^[e^J-i/;. 

Moser's type estimates then imply that for all r > + 1, \£Z\}jr ~ 
and hence |v|jys+i - is bounded from above by C\U\ s+7/2. This is 

also the case of the second component of D^, as a direct consequence 
of (j2.5p . and because l'^/Jl < 1. 

To control the third component of D^, namely, supjQ vTt^ \ b\ijs+i (with 

b given by (|4.5p ). we need to bound \vdtZ\fjs+i from above. Using 
Theorem 13.11 to compute explicitly dtZ_, one finds 

-e^''d,C ■ {eZ)V( - e v7;9iCdiv-,v - ^ffKl (9iC(eZ) 

which is a smooth function of e^/JI < e/v, V^C, dtC, V^dtC, V^dttp, 
eZ_, -^Q[e(]dtip and -^Q[eC](^dtC{£Z_)y The sought after estimate on 
|b|j:^s+i thus follows from Moser's type estimates (note that Remark 
13.21 is used to control --^^[e(^] {dtC{£:Z_)) in terms of Sobolev norms of 

dtC and eZ_). 

The estimate on C_i is obtained exactly in the same way and we omit 
the proof. □ 

□ 



4.2. Proof of Proposition 

Deducing Proposition 14.11 from Proposition 14.31 is only a technical 
step, essentially based on the equivalence of the norms E'^ and \ ■ Ix" 
stemming from Proposition 13. 4i We only give the main steps of the 
proof. 

Step 1. Since U = {Vi, V2 + sZVi), one can expand |f7|x= in terms 
of Vi and V2 and control the different components using the norm E'^ 
to obtain: 

\U\xs <Cx {E^+\V) + {\U\xs+^./2E'''-'\V)),^to)- (4.19) 

Step 2. Using Proposition O to control ^*+i(U) and E*o+'^(y) in 
terms of U^ = {U^,U^ - eZ_\t=oU^) and H = {Gi,G2 - eZGi) in 
(|4.19p . one gets 
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Step 3. Replacing H by (Gi,G2 - eZ_Gi) and by {Ul,U^ - 
eCk=of^i) one obtains the following control on T''^'^{t^V^ ,H) (r = 
s,to): 

T'-+\t, V^H) < C{T+\t, C/O, G) + m ^.+,/,J*«+2(i, u\ G))r>to)- 
Step 4. The proposition follows from Steps 2 and 3. 

5. Main results 

5.1. Large time existence for the water-waves equations 

In this section we prove the main result of this paper, which proves 
the well-posedness of the water-waves equations over large times and 
provides a uniform energy control which will allow us to justify all 
the asymptotic regimes evoked in the introduction. Recall first that 
the semi-normed spaces {X'^, \ • have been defined in Definition 
Olas 

X' := {(C,V'),C e G H'-^^^R^}, 

and \{Ct'4')\x» ■~ ICIh" + l^^l/f^, and define also the mapping a by 

a(C, V) := ^(V^V' - eZ[eC]V\) ■ V^Z[eC]i; 

- eZ[eC]{C + ^2[(C, V')]) + edcZ[e-]i; ■ g[eC]^ + 1, 

where Z[eQ is as defined in Theorem 13. II and A2 is defined in (j4.2p . 
and where = (1 + ^JJI)~^ . The only condition we impose on the pa- 
rameters is that the steepness £^/Ji and the ratio /3/e remain bounded. 
More precisely, (e, 7, (5) G Vm {M > 0) with 



pjvf ={(e, fi, 7, P) G (0, 1] X (0, 00) X (0, 1] X [0, 1] , e^fi < M and ^<M}. 

We can now state the theorem: 

Theorem 5.1. Let to > 1, > and P C Vm. 

There exists P> D >0 such that for all s > sq, 6 G i7*+^(M2), and 
all family (Cp,'0p)pg-p bounded in X^~^^ satisfying 

inf 1+eC- pb>0 and inf a(C°, V'J!) > 

(uniformly with respect to p = (e,^,7, /?) G V), there exist T > 
and a unique family {Cp,'^p)peV bounded in C([0, solving 
(jl.4p with initial condition (Cp,ipp)p^p. 
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Remark 5.1. The time interval of the solution varies with p £ V 
(through £ and v); when we say that (^p, 'i/'p)pe-p is bounded in 
C([0, X'^'^^), we mean that there exists C such that 

yp er, vt G [0, \Cp{t)\H^+D + mp{t)\H^+^ < c. 

Remark 5.2. For the shallow water regime for instance, one has e = 
/3 = 7 = 1 and fi is small (say, ^ < 1); thus, we can take V = 
{1} X (0, 1) X {1} X {!}; for the KP regime (with flat bottom), one 
takes V = {{£,£, ^/£),£ G (0,1)} X {0}, etc. 

Remark 5.3. The condition ini^2 a{Cp,tpp) > is the classical Tay- 
lor sign condition proper to the water-wave equations ( [ISmSlESlIMl 
[55 | [3|[T2 l B5] ■ among others). It is obviously true for small data and 
we give in Proposition 15.11 some simple sufficient conditions on the 
initial data and the bottom parameterization b which ensure that it 
is satisfied. 

Remark 5.4. One also has the following stability property (see Corol- 
lary 1 in |2j): let T > and (C/p^^)pg-p = (Cp^^, ■0p^^)pg-p, bounded in 
y^"*"^ (see Definition 14. 2p . be an approximate solution of (jl.4p in the 
sense that 

with {Rp,rp)p bounded in C([0, X^+^) n C7i([0, if ]; ^ 
X^^^ (and Sp > 0). If moreover the Up^^ are admissible, then one 
has 

yt e [0, ^ inf{T, T}], \Up{t) - U^^P{t)\^^+o < Cst 5p, 

where Up G (^([0, X*"*"^) is the solution furnished by the theorem. 
For 5p small enough, one can take T = T_. 

Remark 5.5. The numbers P and D could be explicited in the above 
theorem (as in Theorem 1 of [2] for instance), but since the focus 
here is not on the regularity of the solutions, we chose to alleviate 
the proof as much as possible. For the same reason, we use a Nash- 
Moser iterative scheme which allows us to deal with all the different 
regimes at once, though it is possible in some cases to push further 
the analysis of the linearized operator and use a standard iterative 
scheme (as shown in [2j| for the shallow- water regime). 

Proof. Let us denote in this proof e = e/z^ and omit the index p for 
the sake of clarity. Rescaling the time hy t t/e and using the same 
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notations as in ()4.ip and (14. 2p . the theorem reduces to proving the 
well-posedness of the IVP 

dtU + + A[U] = 0, 
U\t=o = 

on a time interval [0, T], with T > independent of all the parameters. 
Define first the evolution operator S^{-) associated to the linear part 
of the above IVP. The following lemma shows that the definition 

S'{t)U^ ■.= U{t), with dtU + -CU = {) and U\^^^ = U^ 

(5.1) 

makes sense for all data G X*. 

Lemma 5.1. For ah C/° G X\ S'{-)U° is weh defined in C{[0,T]; X") 
by dSl]). Moreover, for all < t < T, 

Proof. Proceeding as in the proof of Proposition 14.11 (in the very 
simple case U = (0,0)^), one checks that S^{t)U^ makes sense and 
that the estimate of the lemma holds if G X*. 
Now, let us extend this result to data G X*. Let i be a smooth 
function vanishing in a neighborhood of the origin and being constant 
equal to one outside the unit disc, and define, for all 6 > 0, = 
i(\D\/5). The couple {7°'^ := (C°,t^V'°)^ = (C°,<''')^ then belongs 
to X' and U^{t) := S'{t)U^'^ = {(^ (t) , (t))'^ is well defined in X^ 
Since 

\U\t) - U''{t)\ < CiT, -1, \b\^s+v2, ^, ^)\U^'' - U'^''\xs, 

it follows by dominated convergence that {C^)s^o and (W^)5 ^0 are 
Cauchy sequences in C{[0,T]; H^iR"^)). Therefore, (C^) ^ C and 
(*PV'') ^ w in C([0,r];F^(M2)), as 5 goes to 0. Defining V(i) = 
V'" - i /d C{t')dt' and using V'^(i) = iP^'^ - i C\t')dt\ one deduces 
w = qSV', from which one infers that V^- G C([0, T]; F''-V2(k2)2)_ 

/,<5 



Prom the convergence ^p-i/^ — > a; = ^P^/j in i/*(R ) and Proposition 
one deduces also that Q[Q\iIj^ ^ Q[Q]iIj in H^-^I'^{R?). One can thus 
take the limit as 5 ^ in the relation dtC\t) - 7^^[0]V'^(t) = 0, 

thus proving that [Q, ip) G X* solves the IVP (jS.ip . Since the solution 
to this IVP is obviously unique, this shows that S'^{-)U^ makes sense 
in X"" when C/° G X^ 
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The last assertion of the lemma follows by taking the limit when 
(5 ^ in the following expression 

1 p e. 



□ 



We now look for the exact solution under the form U = S''{t)U^ + V, 
which is equivalent to solving 



dtV + i£F + J^[t, V]=h 
V\t=o = {0,Oy 



T (5.2) 



with J^[t, V] := A[S'{t)U'^ + V]- A[S%t)U'^] and h := -A[S%t)U'^]. 
We can now state two important properties satisfied by £ and 
(in the statement below, the notation J^j^^ means that has been 

differentiated i times with respect to time and j with respect to its 
second argument). 

Lemma 5.2. Let T > 0, p = 1 and m = 5/2. Then: 

i. For all s > to, the mapping C : X^'^'^ is well defined and 
continuous; moreover, the family of evolution operators (5'^(-))o<e<eo 
is uniformly bounded in C{[-T, T]- Lin{X'+"' , X')) . 

ii. For all < i < p and 0<i + j<j)-l-2, and for all s > + im, 
one has 

sup |6*.F«[t, U]{Vu . . . , F,-)|s-.™ < Cis,T, |[/|t„+(,+i)^) 
te[o,r] 

j j 

^C^\yk\s+m]\\Vl\to+{i+l)m + \U\s+m I Vfc , 
k=l l^k k=l 

for all U e F'*+™(M2) and (Vi, . . . , Vj) G H'+'^{R'^y. 

Proof, i. The property on S'^{-) follows from Proposition 14.31 with 
U_ = (0,0) (recall that we rescaled the time variable). In order to 
prove the continuity of £, let us write, for all W = (C, '0)"^, 

\cw\x^ < \—g[o]^p\Hs + -\C\h^ + m\H^ 

<\^Q[mHs + c{t)\c\Hs+.. 

One therefore deduces the continuity property on C from the following 
inequality: 
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for ^ > 1, one has the uniform bound ^ < Xj ^fji^ and the inequahty 
is a direct consequence of Proposition 13. 3t for /i < 1, one has i/ ~ 1 
and we rather use Remark 13.31 

ii. Since by definition T{t,\J\ = A[S%t)U^ + U] - A{S%t)U^), it 
fohows from the first point that it suffices to prove the estimates in 
the case i = and with replaced by A. Recall that A is explicitly 
given by ()4.2p and remark that 



1 

^1^ Jo 

= [ -^GiezCKC^Z) + zV^ ■ {CY.)dz, 
Jo VJ^ VJ^ 



where Z and v are as in Theorem 13.11 (with C = C 3.nd if) = ifj). 
The estimates on A are therefore a direct consequence oTProposition 
[331 □ 



The well-posedness of (j5.2p is deduced from the general Nash-Moser 
theorem for singular evolution equations of [2j (Theorem 1'), pro- 
vided that the three assumptions (Assumptions l',2' and 3' in [2j) it 
requires are satisfied. The first two, which concern the linear oper- 
ator C and the nonlinearity J-\t, •], are exactly the results stated in 
Lemma 15.21 The third assumption concerns the linearized operator 
around V_ associated to ()5.2p ; after remarking that 

dt + ^C + dv^[t, ■] = -C(^,^), 

with U_ = ((',■0)"^ = S'^{t)U^ + y_, one can check that this last as- 
sumption is exactly the result stated in Proposition UTTJ provided that 
the following quantity (which is the first iterate of the Nash-Moser 
scheme, see Remark 3.2.2 of [2]) 

Uo:=t^ S'{t)U° + / S'{t - t')T[t\ U^]dt' (5.3) 
Jo 

is an admissible reference state in the sense of Definition 14.11 on the 
time interval [0, T] (recall that we rescaled the time variable). Taking 
a smaller T if necessary, it is sufficient to check the admissibility 
at t = 0, which is equivalent to the two assumptions made in the 
statement of the theorem (after remarking that a{C^,'ip^) = a{t = 0), 
with a as defined in ()4.5p and U_ = Uq). The proof is thus complete. 

□ 

We end this section with a proposition showing that the Taylor 
sign condition 

inf a(Cp, V'p)peP > 0) uniformly with respect to p (5.4) 
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can be replaced in Theorem 15.11 by a much simpler condition. We 
need to introduce first the "anisotropic Hessian" Til associated to 
the bottom parameterization 6, 

^ ._ ( dlb 7292 6 
and the initial velocity potential given by the EVP 



+ 7V5y^° + = 0, -l + (5b<z< eCp, 



Proposition 5.1. Let to > 1, ^ > jind V C Vm', let also b G 
F*o+2(e2)^ (C0,V;°)pe7' be bounded in X^^+^ and (<?°)pep solve the 
EVPs Then, 

i. There exists eo > such that (j5.4p is satisfied if one replaces V by 
'Peo ■= {p = (e, 7, /?) e ez^""^ < eo}; 

ii. If there exist /ii > and 7 G C'((0, 1] x (0, ^ui]) such that for all 
p = (e,^,7,/3) G one has /i < /ii and 7 = 7(6, /i), and if 

-e2/3^W2'(V^° U_i+,J < 1, for all peV, 
then the Taylor sign condition ()5.4p is satisfied. 



Remark 5.6. The first point of the proposition is used to check the 
Taylor condition in deep water regime; in this latter case, one has 
indeed e/v ~ e^/Jl which is the steepness of the wave, the small 
parameter with respect to which asymptotic models are derived. 
The second point of the proposition is essential in the shallow-water 
regime (since e/u does not go to zero as ^ 0). It is important 
to notice that it implies that the Taylor condition is automatically 
satisfied for flat bottoms. 



Remark 5.7. S. Wu proved in |48p49j that the Taylor sign condition 
(15. 4p is automatically satisfied in infinite depth; this result was ex- 
tended in [29] to finite depth with flat bottoms. The result needed 
here is stronger, since we want (j5.4p to be satisfied uniformly with 
respect to the parameters. In the IDH-case, for fiat bottoms, and 
in the particular case of the shallow-water regime, such a result was 
established in [331. 



Proof. As in the proof of Theorem \5.1\ we omit the index p to alleviate 
the notations. 

i. As seen in the proof of Theorem 15.11 one has a{C^,ip^) = a{t = 0), 
where is as defined in (j4.5p (with U_= Uq and Uq given by (|5.3p ). 
Thus, a{C^,ip^) = 1 + and |o|l=o > 1 — eo|.^|L°°- It follows from 
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Lemma l4.ll that for the range of parameters considered here, \b\L°° 
is uniformly bounded on [0,T], so that the result follows when eo is 
small enough. 

ii. Step 1: There exists /io > such that |5.-^p is satisfied for all 
p = (e,/i, 7, /3) S V such that fi < /io- It is indeed a consequence of 
Remark 13.31 that \b\Loc, = 0{yJJi) as /x — > 0; since moreover e = f 
remains bounded, one can conclude as in the first step. 
Step 2. The case fi > fiQ. For all time t, let 'P{t) denote the solution 
of the BVP (|5.5p , with the Dirichlet condition at the surface replaced 
by ^° 1,^,^0 = Mt), where Uo{t) = {Co{t),Mt)) is given by 
Let us also define the "pressure" P as 

_ip := 5,^+ i(£|v^<?|2 + —(9,^)2) + -z. 

Since Uq = (CojV'o) solves (|1.4p at t = 0, one can check as in Propo- 
sition 4.4 of [21] that P{t = d,X,£Q^{X)) = 0. Differentiating this 
relation with respect to X shows that -V^C" • V^P = £\SI^ C^\'^d^P 
on the surface, from which one deduces the identity. 



' 



where dr)P\ „ stands for the outwards conormal derivative associ- 



ated to the elliptic operator ^id^ + j'^fJ.dy + d'^. Expressing (p and its 
derivatives evaluated at the surface in terms of ^, one can then check 
that 

Let us now remark that P solves the BVP 

J ifid'^ + 7^9^ + a2)P = h, -1 + Pb<z< eC°, 

with 

2 ^ V jJLU 

Exactly as in the proof of Proposition 4.15 of [29], one can check that 
h < and use a maximum principle (using the fact that (j5.6p links 
a to the normal derivative of P at the surface) to show that ii g <0 
then there exists a constant c(e, /i, /3) > such that a(to) > c(e, /3)- 
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We thus turn to prove that g < 0. 
Recall that by construction of ^, 

Differentiating this relation with respect to j (j = x,y), one gets 
and using this formula one computes 



-2/ 



since moreover (1 + /?^|V6p)^/^9„(z)|^^_^^^j^ = 1, one gets 5 > if 
the condition given in the statement of the proposition is fulfilled. 
As detailed above, we therefore have a{t = 0) > c(e, /i, /3) > 0; more- 
over, there exists by assumption /ii such that for all p = (e, fi, 7, (3) G 
V, one has /i < ^1; due to the first point of the proposition, Step 1 
and the fact the 7 = j{£,fi), it is sufficient to prove the proposition 
for all the parameters p € "Pi with 

Vi := [eo,l]x[/io,Aii]x7([eo,l]x[/io,/^i])x[0, 1] (eo := (l+v7^)"^eo). 

The dependence of c(e, n, (3) > on e, /i and /3 is continuous and 
therefore, inf[eo^i]x[;,o,Ati]x[o,i] c(e,/^,/3) > 0. □ 



6. Asymptotics for 3D water-waves 

We will now provide a rigorous justification of the main asymptotic 
models used in coastal oceanography. 

Remark 6.1. Throughout this section, we assume the following: 

- P and D are as in the statement of Theorem 15.11 

- stands for the initial velocity potential as in Proposition 15. H 

- The bottom parameterization satisfies b G H''^'^^ (M?); 

- Except for the KP equations, we always consider fully transverse 
waves (7 = 1), but one could easily use the methods set in this paper 
to derive and justify weakly transverse models in the other regimes. 
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6.1. Shallow-water and Serre regimes 

We recall that the so-called "shallow-water" regime corresponds to 
the conditions /i <C 1 (so that ~ 1) and e = 7 = 1; we also consider 
bottom variations which can be of large amplitude (/? = 1). With- 
out restriction, we can assume that = 1 (which corresponds to the 
nondimensionalization (jl.4p ). The shallow- water model - which goes 
back to Airy [1] and Friedrichs [19]- consists of neglecting the 0(/i) 
terms in the water-waves equations, while the Green-Naghdi equa- 
tions [21 1 122 ^ 146] is a more precise approximation, which neglects only 
the 0(/i^) quantities. The Serre equations [3ilH6] are quite similar to 
the Green-Naghdi equations, but assume that the bottom and surface 
variations are of medium amplitude: e = (3 = yjjl. 

6.1.1. The shallow-water equations The shallow water equations are 

ra,i^ + vc + iv|y|2 = o, 

\atC + v-((i + c-6)^)=o, ^^-'^ 

and the following theorem shows that they provide a good approxi- 
mation to the exact solution of the water-waves equations. 

Theorem 6.1 (Shallow- water equations). Let s > t^ > 1 and 

(C|J) V')^)o<^t<i be bounded in X'^^^ . Assume moreover that there exist 
/lo > and /^o > such that for all ^ E (0, /xq), 

inf(l + Ql-h)>h, and - /.W^(V^° |^_,^J < 1. 
Then there exists T > and: 

1. a unique family {Cfj.,'4'fj.)o<fj.<fio bounded in C{[0,T]; X''^^^) and 
solving (jl.4p with initial conditions {C'ji^^^)o<fl<^ioy 

2. a unique family (V^^, C^'^)o<M</to bounded in 
C([0,r];i?*+-P-^/2(i^2)3^) j^^j solving I^J^ with initial conditions 

(Cj,VV'°)o<;.<Mo- 
Moreover, one has, for some C > 0, 

VO < ^ < ^0, IC/i-C^^lL°°([0,T]xR2) + IW/i-^/f^lL°°([0,T]xR2) < Cfl. 

Remark 6.2. The existence time provided by Theorem 15.11 is 0(1), 
but is large in the sense that it does not shrink to zero when ^ — > 0. 

Remark 6.3. Instead of assuming that the initial data {C^, V')J)o</i<i 
are bounded in X^^^, we could assume that {C^,^'ip^)o<fj.<i is bounded 
in iJ*+^(IR^)^ (because {^ipl^s+p < {Vipljjs+p, uniformly with re- 
spect to G (0, 1)). 
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Remark 6.4. The 2DH shallow-water model has been justified rig- 
orously by Iguchi in a recent work [24], but under two restrictions: a) 
The velocity potential is assumed to have Sobolev regularity which 
implies that the velocity must satisfy some restrictive zero mass as- 
sumptions and b) the theorem holds only for very small values of /x. 
These assumptions are removed in the above result. 

Proof. The assumptions allow us to use Theorem 15.11 and Proposition 
IS.ll with V = {1} X (0,/xo) X {1} X {1}, which proves the first part of 
the theorem. 

The second point of the theorem is straightforward since 
\Vip^\j^s+p-i/2 < l^fp^lno+p (recall that /u < 1 here), and because 
(j6.ip is a quasilinear hyperbolic system (since inf]g2(l + ( — b) > 0). 
In order to prove the error estimate, plug the expansion furnished by 
Proposition l3.8l into ()1.4p and take the gradient of the second equation 
in order to obtain a system of equations on C^t and = Vtp^. One 
gets 

dtC^, + V-{{l + C^.-b)V^)=fiRl, ^^-^^ 

with {Rl,Rl) uniformly bounded in L°°([0, T]; i?*o(]R2)2+i). An en- 
ergy estimate on (|6.2p thus gives a Sobolev error estimate from which 
one deduces the L°° estimate of the theorem using the classical con- 
tinuous embedding i/*" C □ 

6.1.2. The Green-Naghdi and Serre equations Though correspond- 
ing to two different physical regimes, the Green-Naghdi and Serre 
equations can both be written at the same time if one assumes that 
e = 1 for the Green-Naghdi equations and e = y^I for the Serre 
equations in the formulation below: 

{h + ^T[h, eb])dtV + hVC + ■ 



lV{h^Vvdw{V)) + Q[h,eb]{V) =0 (6.3) 



dtC + v- (hv) = 0, 

where h := 1 + e(C — b) while the linear operators T[h, b] and T>v and 
the quadratic form Q[h,b]{-) are defined as 

T[h, b]V := -^V{h^V ■ V) + ^ [V{h^Vb ■ V) - h^VbV ■ V] 

+hVbVb • V, 
Vv := -(y • V) +div(y), 

Q[h, b]{V) := ^V{h^iV ■ Vfb) + h{^Vvdiv{V) + (V ■ V)%)Vb. 

Both the Green-Naghdi and Serre models are rigorously justified in 
the theorem below: 
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Theorem 6.2 (Green-Naghdi and Serre equations). Let s > 

to > 1 and (C^) V')^)o</i<i be bounded in X*"*"^. Let e = 1 (Green- 
Naghdi) or e = y/]2 (Serre) and assume that for some ho > 0, fiQ > 
and for ah /x G (0, //q), 

inf(l + e(C° - b)) > ho and - f^e^H^iVPl u_,^,J < L 

Then there exists T > and: 

1. a unique family (C/^; V'At)o<^t<^to 

bounded in C([0, f ]; X"+^) and 
solving (jl.4p with initial conditions {C'ji^''P^)o<fi<fio'j 

2. a unique family (y^^^, C^^)o</.<mo bounded in C7([0, //^(R2)3) 
and solving (j6.3p with initial conditions (C^, (1 — e6])VV'JJ) 
(with /i° = 1 + e(C° - &)). 

Moreover, one has for some C > independent of ^ G (0, /io); 

2 

lCM-C^^lLoo([o,i:]xM2)+|VV'/.-(l+^'^[^,e&])V'^^^lioo([o,z:]xR2) < C'y- 

Remark 6.5. The precision of the GN approximation (e = /? = 1) 
is therefore one order better than the shallow-water equations. This 
model had been justified in IDH and for flat bottoms by Y. A. Li 
[33]. The theorem above is stated in 2DH but one can cover the open 
case of IDH non-flat bottoms with a straightforward adaptation. 

Remark 6.6. In the Serre scaling, one has e = ^/JI, and the precision 

of the theorem is therefore 0(/i^/^), which is worse than the 0{iJ?) 
precision of the GN model, but the approximation remains valid over 
a larger time scale (namely, 0(;U~^/^) versus 0(1) for GN). Notice 
also that at first order in n, the Serre equations reduce to a simple 
wave equation (speed ±1) on C and V, which is not the case for GN 
where the shallow-water equations (|6.ip are found at first order. 

Proof. The first assertion of the theorem is exactly the same as in 
Theorem 16.11 in the GN case. For the Serre equations, it is also a 
direct consequence of Theorem 15.11 and Proposition 15.11 with V = 
{(\/i"> ^, 1, \/^), e (0, fio)}- 

For the second assertion, we replace Gfj.,'y[£:Ci Pb] in (jl.4l) by the expan- 
sion given in Proposition 13.81 and take the gradient of the equation 
on ^ to obtain 

/ dtVf, + VC^ + |V|y^|2 - f V(/iV • - eVb ■ V^)^ = fi^R},, 
\dtC^ + V-{hV)=fi^Rl, 

(6.4) 

with {Rl,Rf^)^ bounded in L°°([0, |]; i7*o(]R2)2+i) ^^ile V is defined 
as V := V^- i^T[h,b]V^, so that = V + ^T[h,b]V + 0(//2). 
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Replacing by this expression in ()6.4p and neglecting the 0(/i^) 
terms then gives (16.3p . The theorem is then a direct consequence of 
the well-posedness theorem for the Green-Naghdi and Serre equations 
proved in [2] and of the error estimates given in Theorem 3 of that 
reference. □ 



6.2. Long-waves regime: the Boussinesq approximation 

The long- wave regime is characterized by the scaling 7 = 1, /i = e <C 
1, so that one has ~ 1. As for the shallow- water equations, we 
take u = \ for notational convenience. When the bottom is non-flat, 
it is assumed that its variations are of the order of the size of the 
waves, that is, (3 = e. Since the pioneer work of Boussinesq [8], many 
formally equivalent systems, generically called Boussinesq systems, 
have been derived to model the dynamics of the waves under this 
scaling. Following [7j, these systems where derived in a systematic 
way in [61151110^19]. In [5l[9] some interesting symmetric systems where 
introduced: 

' (1 - ea2A)dtV + VC + e{^V\V\^ + ]^{V ■ V)V + ^FV • V 

Se,puP2 I +^V|C|' - \bVC + aiAVC) = 0, 

^ (1 - ea4A)dtC + V-y + |(V- ((C - b)V) + a-^AV ■V) = 0, 

where the coefficients aj {j = 1,...,4) depend on pi,p2 G M and 
9 G [0, 1] through the relations oi = — ^2 = — Pi); 

as = ^^^-p2, and 04 = ^-^(1 — P2)', some choices of parameters 
yield oi = 03 and 02 > 0, 04 > 0, and the corresponding systems 
S'q p2 are the completely symmetric systems mentioned above. The 
so-called Boussinesq approximation associated to a family of initial 
data (C°,V'°)o<e<i is given by 

C^^ = C^ and F/^'^' = (l-£(l-02)^)(l-£(Cf-6)y/), (6.5) 

where {V^^ ,C^)o<e<i solves Sg^^^^^ with initial data 

y/.°=(l + |(C°-6))(l-|(l-02)Z^)-Vv^o and Cf '° = C°. 

(6.6) 

The following theorem fully justifies this approximation. 

Theorem 6.3 (Boussinesq systems). Let s >tQ > 1 and 

{Ce, V'e )o<e<i be bounded in X'"^^^ and assume that there exist Hq > 
and Eq > such that for all e G (0, Eq), 

inf(l + e(C° - b)) > ho and - e'U^iV^l |^_,^^J < 1. 
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Then there exists T > and: 

1. a unique family (Ce; '0e)o<£<£o bounded in C{[0,j];X'^~^^) and 
solving (II. 4p with initial conditions (Ce > '0e)o<e<eo; 

2. a unique family {V^^ ,Cf)o<e<eo bounded in C([0,^]; F'^+^'s (M^)^) 
and solving S'q with initial conditions (j6.6p . 

Moreover, for some C > independent of e G (0,eo), one has 

VO < t < |, lCe(t) - C^P(t)|oo + |VVe(i) - v;"^p(t)|oo < Ce\ 
where (V^"^^,Ce^^) is given by (j6.5p . 

Remark 6.7. The above theorem justifies a// the Boussinesq systems 
and not only the completely symmetric Boussinesq systems consid- 
ered here: it is proved in [5,9j that the justification of all the Boussi- 
nesq systems follows directly from the justification of one of them, in 
the sense that their solutions (if they exist!) provide an approxima- 
tion of order 0{e^t) to the water-waves equations. 

Proof. The "if-theorems" of prove the result assuming that the 
first statement of the theorem holds, which is a direct consequence of 
Theorem 15.11 and Proposition 15.11 with V = {(e,e, l,e),e G (0,eo)}. 

□ 

6.3. Weakly transverse long-waves: the KP approximation 

We recall that the KP regime is the same as the long-waves regime, 
but with 7 = y/e. Moreover, we assume here that the bottom is 
flat, for the sake of simplicity. The KP approximation [25] consists in 
replacing the exact water elevation by the sum of two counter prop- 
agating waves, slowly modulated by a KP equation; more precisely, 
one defines C,^^ as 

Cf^(t,x) = l(C+(rf,V^y,2;-t) + C-(ei, V^y,x + i)) (6.7) 
where C,±{t,Y,X) solve the KP equation 

This approximation is rigorously justified in the theorem below: 

Theorem 6.4 (KP equation). Let s > to > 1 and (C°,'0°) S 

and assume that there exist /iq > and eo > such that for all 

£ e (0,eo), 

inf(l + eC°) > ho, 

and assume also that (S^S^V", ^^C") G d'^H'+P {R'^f . 
Then there exists T > and: 
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1. a unique family (Cei '0e)o<e<eo solving (II. 4|) with initial conditions 

(C°,'0°) and such that {Ce)o<e<eo, {dx'4>s)o<s<so and {^/£^y■^pe)o<e<eo 
are bounded in C([0, H'+^-^/^); 

2. a unique solution C± G C{[0,T]; H'+p-^/^{R'^)) to (KP)± with 
initial condition {(^ it dxil^^)/2. 

Moreover, one has the following error estimate for the approximation 



Remark 6.8. The very restrictive "zero mass" assumptions that 
dydxip^ and dy(^ are twice the derivative of a Sobolev function comes 

from the singular component d'^dy of the KP equations (KP)-i-. Fur- 
thermore, the error estimate is much worse than for the Boussinesq 
approximations. These two drawbacks are removed if one replaces 
the KP approximation by the approximation furnished by the weakly 
transverse Boussinesq systems introduced in [32j. As shown in [32j, 
the first assertion of the above theorem rigorously justify these sys- 
tems: they provide an approximation of order O(e^t) on the time 
interval [0, T/e], and do not require the "zero mass" assumptions. 

Proof. As for the Boussinesq systems, we only have to prove the first 
assertion of the theorem, and the whole result then follows from the 
"if-theorem" of 02j. Taking V = {(e, e, ^e, 0), e G (0,eo)}, Theo- 
rem [5TT] and Proposition 15.11 give a family of solutions (Ce; '0£)o<e<eo 
bounded in C([0, |]; In particular, (l^pV'el H!>+D)e is bounded, 

and thus (|V''''i/'e|/^s+D-i/2)e is also bounded. Since 7 = -y/e, one has 
\dx'4>\H!>+D-i/2 + \/e\dy'ilJe\H'+D-i/2 < I V'^'i/'e|_H-s+D-i/2 and the claim 
follows. □ 



6.4- Deep water 

6.4-1- Full dispersion model We present here the so-called full dis- 
persion (or Matsuno) model for deep water-waves. Contrary to all 
the asymptotic models seen above, the shallowness parameter ^ is 
allowed to take large values (deep water) provided that the steepness 
of the waves e^/Jl remains small; without restriction, we can therefore 
take V = here (i.e., we use the nondimensionalization ()A.2p ). 

Introducing e = the full dispersion model derived in [36 | [37 1 [TT] 
can be written in the case of flat bottoms (/3 = 0): 



(16^1): 



lim|C.-C 



-KP\ 

lL°°([0,f ]xIR2) 



= 0. 



dtC - T^V + e(T^(CVT^y) + V • (c^)) = 0, 
dtV + vc + e(iv|y|2 - VCT^VC) = 0, 



(6.8) 
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where 7^ is a Fourier multiplier defined as 

Since /? = (flat bottom) and 7=1 (fully transverse), the full 
dispersion model depends on two parameters (e, fi) which are linked 
by a small steepness assumption: 

3eo > 0, (e, fi) e V^o C {(e, /i) e (0, 1] x [1, cx)), e := < eo}. 

The well-posedness of the full-dispersion model has not been investi- 
gated yet, but we can prove that if a solution exists on [0, =] (e > 
small enough), then the solution of the water-waves equations exists 
over the same time interval and is well approximated by the solution 
of the full-dispersion model: 

Theorem 6.5 (Full-dispersion model). Let eo > 0, Q > P large 
enough, s > to > 1 and € X''^^^ , and assume that 

VeG(0,l], inf(l+e(C°-6)) > /lo > 0. 

Let also T > and let {C^^j^, K^^)(e,^)ep,Q be bounded in 

C([0, =],H'+'^{R'^f) and solving (USD with initial condition (C°, V^°- 

e(T^V>)VC°). 

Then, if eo is small enough, there is a unique family (Ce,/^, '0e,At)(e,/x)e'PEQ 
bounded in C([0, =];X^~^^) and solving p.4|) with initial conditions 
In addition, for some C > independent of {e,fi) G V^q, one 

has 

ICs./.-cil'/f li«=([o £]xr2)+IW^-'^/^Il-([o £]xM2) ^ •= ^\/^)- 

Proof. Since 7^ : i7''(M^)^ H''"{M?) is continuous with operator 
norm bounded from above by 1, the mapping V G i?''(M^)^ t-^ V — 
e{T^V)VC G i/^(M2)2 is continuous for ah r > to and C G H''+\R'^). 
Moreover, this mapping is invertible for e small enough, and one can 

accordingly define V := {I- eVCT^)-^V , so that V = V-e(Tf,V)VC- 
Replacing V by this expression in (|6.8p gives 

f atC - T^v + 6(v • icv) + r^vicr.v)) = eV^, 

\ + VC + ei(V|y|2 - V(T^y)2) = e^Vr^, 

and where the exact expression of := (r^,Vr^)e is of no impor- 
tance. Now, let dt + L denote the linear part of the above system and 

S{t) its evolution operator: L := ^ ~ (C)^)) 
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S{t)U := u{t), where u solves {dt + L)u = 0, with initial condition 
u\t^o = U. Since 7^ is a Fourier multiplier, one can find an explicit 
expression for S{t), but we only need the following property: S{t) is 
unitary on (r S R) defined as 

■={U = {C,V) G H^iR^f, 

\TT\ lAi , |/tanh(^|^|) 1/2 I 

\U\zr := \C\h- + K ^ ) ^1^,. < ooj. 

Writing u := ((", V), we define w := {(, W) as 

w-u-e^ [ S{t-t')R,{t')dt'; 
Jo 

remarking that \V\Hr-i/2 < | ( '^"^j^'^' Y ^V|^, and that 
yf^m'+Vro2^ . tanh(^|g|) 1/2^ 

uniformly with respect to > 1, and since S{t) is unitary on Z^, one 
gets 

Vr > 0, sup \w{t) - u{t)\H^ < eT sup(|r,^|j:^,+i/2 + |r^|jyr+i). 
[0=] [o,f] 

(6.10) 

Furthermore, one immediately checks that w solves 
jdtC-T^W + ef\C,W) = e^kl 

2/A - .2y71,2 



dtW + VC + eVf\C, W) = e^Vk, 



with initial condition w\^_^ = {(^, Vt/j^)'^, and where := ^{f^{w) — 
r{u)), := \{f\w) ->(S)), and 

f\C,V) := V • (CF) +T^V(CT^y), /2(C,y) := \{\V\^ - {T^Vf)- 

from (j6.10p . one gets in particular that 
\{klkl)\^s^P + \{dtkldtkl)\^.^+P.,/, <C{T^{^^^^^ (6.11) 

provided that Q is large enough. 

Using the fact that all the terms in the equation on W are a gradient 
of a scalar expression, as well as M^i^^q, it is possible to write w = 

(C,VV')'^, and (CV') solves 

f dtl-T^Vi^ + e/HC, VV) = e2A;,\ 



(6.12) 

5tV' + C + e/'(C,VV) = e2fe2, 
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with initial condition {C,tp)\^^Q = {C^,tp^). 

Remarking now that ^[O]-;/' = ^/JlT^Vip and writing U = (CjV')) oris 
can check that (16.121) can be written 

dtU + CU + eA^^^[U] = e\kl, , (6.13) 

where A^^^ [U] is given by the same formula (j4.2p as ^[C/], but with 
the Dirichlet-Neumann operator QleQip replaced by the first order ex- 
pansion given in Proposition l3.9l and with the O(e^) terms neglected. 
One thus gets 

dtU + CU + eAp] = e^H„ 

withFe = {kl,k'^f+l{A^^^[U]-A[U]). From (IHTH . PropositionESl 

and (|6.1ip . one has (-ffe)e is uniformly bounded in C([0, n 

we can therefore conclude with Remark 15.41 
and Proposition 15.11 □ 

6.4-2. A remark on a model used for numerical computations The 
Dirichlet-Neumann operator is one of the main difficulties in the nu- 
merical computation of solutions to the water-waves equations (jl.4p 
because it requires to solve a d-l- 1 (d = 1, 2 is the surface dimension) 
Laplace equation on a domain which changes at each time step. A 
common strategy is to replace the full Dirichlet-Neumann operator 
by an approximation which requires less computations. An efficient 
method, set forth in [15j, consists in replacing the Dirichlet-Neumann 
operator by its n-th order expansion with respect to the surface ele- 
vation C,. When n = 1, it turns out that the model thus obtained is 
exactly the same as the system (I6.13P used in the proof of Theorem 

EH 

We can therefore use this theorem to state that: the precision of 
the modelization in the numerical computations of [iS] is of the same 
order as the steepness of the wave. 

One will easily check that when the n-th order expansion is used, then 
the precision is of the same order as the n-th power of the steepness. 



A. Nondimensionalization(s) of the equations 

Depending on the value of fi, two distinct nondimensionalizations are 
commonly used in oceanography (see for instance [iBj). Namely, with 
dimensionless quantities denoted with a prime: 

— Shallow- water, ie /i <C 1, one writes 

X = Ax', y = -y' z = dz' t = -^t' , . , 

Q = aC, <P = ^XVgd^', b = Bb'. ^ ' 
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— Deep-water, ie fi ^ 1, one writes 

X = Xx', y = ^y', z = Xz' , t = -^t' , 

( = aC, <P = ay/^<^', b = Bb'. ^ ' ' 

Remark that when ~ 1, that is when A ~ both nondimensional- 
izations are equivalent, we introduce the following general nondimen- 
sionalization, which is valid for all /i > 0: 

X = Ax', y = ^y', z = dvz', t = -^^t' , 



C = aC', ^=!iX^JfcP',b = Bb', 

where v is a smooth function of /u such that ~ 1 when ;U ^ 1 and 
u ~ ^"^/^(= X/d) when /i > 1 (say, u = {I + ^)~^). 

The equations of motion (II. ip then become (after dropping the 
primes for the sake of clarity) : 

+ 1/^7 V^^^ + = 0, -{-1 + Pb) < z < -c, 
-u^fiV^i^b) ■ V^<P + d,<P = 0, z = i(-l + (3b), 
9tC - A ( - ^Vv^(-C) • V^^ + d^cp) =0, z = -c, 

Q,^ + ^(^iv-^p + -^,{m') + c = 0, z = £c. 

(A.3) 

In order to reduce this set of equations into a system of two evo- 
lution equations, define the Dirichlet-Neumann operator ^^^^[^Cj 
as 



with ^ solving the boundary value problem 

v^lidl^ + v'^-i'^[idl<^ + dl<l> = 0, -{-I + pb)<z< -C, 

(as always in this paper, stands for the upwards conormal deriva- 
tive associate d to the elliptic equation). As remarked in [52 t [T7 l [T6] . 
the equations ()A.3P are equivalent to a set of two equations on the free 
surface parameterization ^ and the trace of the velocity potential at 
the surface V = ^1^-^/^^ involving the Dirichlet-Neumann operator. 
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Namely, 



5,^ + C + ^|v-^l 2(i + .V|v^cP) = 

(A.4) 

In order to derive the system p.4p . let ^^^^[eC? be the Dirichlet- 
Neumann operator ^^^^[^Ci Pb]- corresponding to the case v = 1. One 
will easily check that 

1/ ' 1/ 

so that plugging this relation into ()A.4p yields 

a,^ + C + ^|V-^| 2(l + eV|VXP) = 
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